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ABSTRACT 



Keywords: Wave packets; Kerr medium; revivals; fractional revivals; co- 
herent states; photon-added coherent states; expectation val- 
ues; higher moments; squeezing; Wigner function; entropy; en- 
tanglement; time series; Lyapunov exponent. 

Starting with an initial wave packet that corresponds to an ideal coherent 
state, we show that, in the framework of a generic Hamiltonian that models the 
propagation of a single-mode electromagnetic field in a Kerr-like medium, distinc- 
tive signatures of wave packet revivals and fractional revivals are displayed by the 
time evolution of the expectation values and cumulants of appropriate observables. 
This enables selective identification of different fractional revivals. Using a class 
of photon-added coherent states as the initial states, we further show that the 
revival phenomena provide clear signatures of the extent of the departure from 
coherence, and of the deviation from Poissonian photon number statistics, of the 
initial state of the radiation field. A comparison is also made of the manner in 
which other non-classical effects are displayed as the wave packet evolves, in the 
cases of an initial photon-added coherent state and an initial ideal coherent state, 
respectively. These effects comprise squeezing, higher-order squeezing, and the 
non-positivity of the Wigner function. 

Extending these studies to the case of two interacting modes, we consider 
another generic Hamiltonian that models the interaction of a single-mode electro- 
magnetic field propagating in a nonlinear medium, an atom of which is modeled 
by an anharmonic oscillator. The dynamics of expectation values and higher mo- 
ments of observables pertaining to the field sub-system as well as the total system 
are studied over a range of values of a relevant parameter, namely, the ratio of 
the strength of the anharmonicity to that of the inter-mode coupling. Analogs of 
revival phenomena are found, and their signatures identified in quantifiers such 
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as sub-system entropies that serve to characterize the degree of entanglement. A 
time series analysis of the mean photon number (or, equivalently, the mean energy 
of the field mode) is carried out, to deduce the embedding dimension of the recon- 
structed phase space and the maximal Lyapunov exponent. The corresponding 
dynamical behavior is found to range from mere ergodicity (as in quasi-periodicity) 
to exponential instability characterized by a positive maximal Lyapunov exponent, 
as the parameter ratio mentioned above is varied. 

Our results shed light on the non-classical effects that occur during the prop- 
agation of a quantum mechanical wave packet in a nonlinear medium. Several 
avenues for further research are also revealed. 
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CHAPTER 1 



Introduction 



The long-time behavior of classical dynamical systems, and the degree of random- 
ness they exhibit (ranging from mere ergodicity to fully-developed chaos), have 
been investigated extensively in the literature. They continue to be subjects of 
abiding interest and current research. Ergodicity theorems of wide applicability 
ensure that, generically, the dynamics is metrically transitive — either on the 
energy surface, in conservative (Hamiltonian) systems, or on some attractor, in 
dissipative systems. While the evolution equations for the dynamical variables 
describing the behavior of a classically chaotic system are necessarily nonlinear, 
the corresponding quantum dynamics is governed by the Schrodinger equation 
(or, in the more general case of mixed states, by the Liouville equation), which 
is inherently linear. It is not straightforward, therefore, to assess the extent of 
randomness in a quantum system, or the precise manner in which the information 
content of a quantum system changes with time. 

At least two different approaches have been employed in identifying signatures 
of the nature of the ergodicity of the quantum mechanical counterpart of a generic 
classical system [p]]. One approach relies on the observation that such signatures 
are generally manifested in the statistics of the energy levels of the corresponding 
quantum system. If the system is classically integrable, the quantum levels cluster 
together, and could even cross when a specific parameter in the Hamiltonian is 
varied ||. A classically chaotic system, on the other hand, has its corresponding 
quantum energy levels so correlated as to resist such crossings 0, ffl. Another 
approach is based on investigating the dynamics of the overlap between two quan- 
tum states of the same physical system which originate from the same initial state, 
but with slightly different values of one of the control parameters 0. The time- 
dependent overlap is close to unity at all times if the normalized initial state is 
located in a classically regular region of phase space. In contrast, if the initial state 



is in a chaotic region of the classical phase space, the overlap falls off exponentially 
in time. 

The foregoing lines of investigation concern quantum signatures of classical 
dynamics. The inverse problem, namely, the identification of signatures of non- 
classical effects in the temporal behavior of quantum mechanical expectation val- 
ues (which, in turn, could be regarded as effective classical dynamical variables in 
an appropriate phase space), has also received attention. The dynamics of a quan- 
tum wave packet governed by a nonlinear Hamiltonian provides adequate scope 
for such an investigation to be carried out, as a wide variety of non-classical effects 
such as revivals and fractional revivals (for reviews see, e. g., 0, [7|, §), as well 
as squeezing (for a review see, e. g., ||), are displayed by the wave packet as it 
evolves in time. 

While a generic initial wave packet (corresponding to the state |^(0))) governed 
by a nonlinear Hamiltonian spreads rapidly during its evolution, it could return to 
its original state (apart from an overall phase) at integer multiples of a revival time 
T rev , under certain conditions. This is signaled by the return of | (ip(0)\ip(t)) | 2 to 
its initial value of unity at t = nT rev . The role played by the non-quadratic terms in 
the quantum Hamiltonian on the dynamics of a wave packet as it spreads and loses 
its original shape is best understood by examining the dynamics of appropriate 
expectation values and the consequences of Ehrenfest's theorem. Consider the 
simple one-dimensional example of a particle of mass m subject to a potential 
V(x). The Hamiltonian is 



where x is the position operator and p is the momentum operator. Ehrenfest's 
theorem states that 





d (x) (p) 




dt m 



dip) 



dV{x) 




dt 



dx 
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For quadratic Hamiltonians, the latter equation gets simplified to 

dip) __ dV{(x)) 
dt dx 



;i.4) 



Correspondingly, these expectation values follow the classical motion. For non- 
quadratic Hamiltonians, (approximate) classical motion arises only if the potential 
can be approximated to the quadratic form. In fact, certain techniques for the 
semi-classical treatment of the propagation of wave packets HlO| , |TT| are based 
on the replacement of the exact Hamiltonian with a tractable, usually quadratic, 
approximation governing the evolution of the wave packet. Further, the quadratic 
potential turns out to be the most general potential under which a Gaussian wave 
packet continues to remain a Gaussian. In fact, there is an interesting parallel 
jl2fl between the following two situations: The quadratic approximation given in 
Eq. (|1.4j) , and the retention of the Gaussian shape under a quadratic potential, on 
the one hand; and the role of the second-order term in a series expansion for the 
classical Hamiltonian in the linearized behavior of classical orbits close to some 
reference orbit in phase space, on the other. 

The eventual spreading of a generic wave packet (which could initially follow 
classical motion for a brief period of time) is inevitable if a quadratic approxima- 
tion does not hold good. It has been argued |T^| that such wave packet spreading 



in both position and momentum space should also be describable in a classi- 
cal framework by considering an ensemble of nearby classical orbits. Despite the 
spreading phenomenon, interesting quantum interference phenomena [13[ can lead 
not only to wave packet revivals, but also to fractional revivals between two suc- 
cessive revivals at specific instants. The latter is characterized by the splitting 
up of the initial wave packet into a number of spatially distributed sub-packets, 
each of which is similar to the original wave packet. Classical analogs of revival 
phenomena can also be envisaged. For instance, one can model the individual 
energy eigenstates comprising the wave packet as an ensemble of racing cars on a 
circular track with different speeds |14j] . The quantum mechanical spreading of the 
wave packet arises due to the difference in speeds. However, certain distinct time 
scales appear in the problem: The initial classical periodicity may be retained for 
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a few laps, while for longer times no correlations between the various racers are 
observed — but obvious patterns can recur, including the original form (revivals), 
as well as small similar groups of racers clumped together at certain instants (frac- 
tional revivals). Another interesting analogy exists between wave packet revivals 
(in time) and the recurrence of images (in space) in the classical Talbot effect. 
This is the repeated self-imaging of a diffraction grating at multiples of a certain 
fundamental distance from the source due to coherent interference of waves, akin 
to quantum revivals of a wave packet at multiples of the fundamental time T rev . 
Further, in between the appearance of two such self-images, Talbot images (which 
are superposed copies of the original) appear at definite spatial positions, akin to 
fractional revivals of the quantum wave packet at certain instants between two 
successive revivals (see, for instance, |L5|, |l6ll). 

In this thesis, we focus on the quantum revival phenomena and other non- 
classical effects that arise close to instants of revivals and fractional revivals, and 
the information that can be obtained about these effects by examining the expec- 
tation values and higher moments (or cumulants) of operators corresponding to 
appropriately chosen observables. It is evident that, at full revivals of any initial 
state, the expectation values of all observables return to their initial values. Thus, 
in a "phase space" where these expectation values play the role of classical dynam- 
ical variables, the trajectories are closed orbits, indicative of the regular temporal 
evolution of the corresponding dynamical system. In principle, such a classical 
phase space would be infinite-dimensional, because the full information contained 
in a specific quantum state can only be obtained if the expectation values and all 
the higher moments of all operators pertaining to the system are included in the 
set of dynamical variables. A more practicable approach, however, is to use the 
time series of an available set of expectation values pertaining to the system, in 
order to estimate the embedding dimension of the relevant phase space, as is done 
when carrying out a time-series analysis to determine the extent of randomness of 
a given set of signals. 

Interesting signatures of different fractional revivals also manifest themselves in 
the dynamics of the expectation values of various operators. For one- dimensional 
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systems with Hamiltonians that can be written in the form given in Eq . (|1 . 1| ) , the 
Ehrenfest equations can be generalized to include the temporal evolution of the 
second moments of x and p. These are given by 

d (x 2 ) 1 

__ = _ ( X p + px) , (1.5) 



= *M = Ml + {xF{x)) , (i. 6) 
dt dt m 

and 

^l = (pF(x) + F(x)p), (1.7) 

where F(x) = —dV/dx. These equations have been used to examine the role of 
quantal uncertainties in wave packet dynamics in the case of the free particle and 
the harmonic oscillator \T%. We have examined the role of expectation values in 



a more general setting by considering a nonlinear Hamiltonian which is not of the 
form given in Eq. ( |1.1|) . We have established that, by tracking the time evolution of 
various moments of certain operators, selective identification of different fractional 



revivals can be achieved 18 



While revival phenomena have been investigated extensively in the context of 



Rydberg atoms, both theoretically and experimentally (see, e. g., [19], [21], gj), 
for our purposes we have studied the dynamics of an initial wave packet (a single 
mode of an electromagnetic field) propagating in a Kerr-like medium. This is 
motivated by the fact that non-classical effects displayed by optical fields have 
also been investigated intensively for well over a decade now. The model that we 
consider (in the first part of this thesis) neatly captures the important role played 
by expectation values of higher moments of operators in wave packet dynamics: 
we examine the dynamics of the field mode under an effective Hamiltonian 

H = h X a^ 2 a 2 , (1.8) 

where a and are photon destruction and creation operators. The starting point 
in obtaining this effective Hamiltonian [23|, which is written in terms of the field 
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operators alone (without explicitly invoking the field-medium interaction terms), 
is the set of classical field equations pertaining to the optical Kerr effect, namely, 
the dependence of the refractive index of the medium on the intensity of the 
light propagating through it. While this effect, displayed by certain nonlinear 
optical materials (Kerr media), can be explained within a classical framework, 
the field amplitudes are replaced in the standard manner by photon destruction 
and creation operators, in order to take into account the quantum nature of the 
radiation field. The effective quantum Hamiltonian is constructed subject to the 
requirement that the Heisenberg equations for the field operators yield the correct 
dynamical equations in the classical limit. A 'Kerr-like' medium is one for which 
the above effective Hamiltonian is a good approximate Hamiltonian. In this thesis, 
we shall use the words 'Kerr medium' and 'Kerr-like medium' interchangeably. 

We note, in passing, that fractional revivals of an initial Gaussian wave packet 
propagating under the Kerr Hamiltonian have been considered to be potentially 
useful in quantum computation. They play a significant role in the context of 
quantum cloning with continuous variables |24| . Fractional revivals corresponding 



to the appearance of two and four superposed sub-packets are particularly useful 
in implementing one- and two-bit logic gates f25|. Further, by combining these 
superposed states with the vacuum using a 50% beam splitter, the output state 



can be shown to be an entangled state (a Bell state) under certain conditions p6 



Another aspect which we have investigated in some detail in this thesis, in 
the context of wave packet evolution in the Kerr-like medium, is the precise role 
played by the initial state in determining the non-classical features exhibited by 
the time-evolved state. Taking the standard oscillator coherent state (abbreviated 
as CS throughout this thesis) as the reference initial state for the purposes of 
comparison, we have examined the extent to which the departure from coherence 
of different initial states affects the dynamics of the wave packet. In particular, 
we have identified, both qualitatively and (wherever possible) quantitatively, the 
manner in which effects such as revivals, fractional revivals and squeezing of the 
state at certain instants depend on the degree of coherence of the initial state. The 
initial states considered for this purpose are the photon-added coherent states (ab- 
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breviated as PACS throughout this thesis). By repeated application of the photon 
creation operator on a CS, photons can be added to it systematically, and a quan- 
tifiable departure of the state from perfect coherence obtained as a consequence. 
The recent experimental production and complete characterization of the single- 
photon added coherent state by quantum state tomography |27J adds impetus 
to our investigation. The single-photon added coherent state is particularly in- 
teresting, as it exhibits marginal departure from coherence, and is intermediate 
between a one-photon Fock state (that is purely quantum mechanical) and a CS 
(a minimum uncertainty, "classical" state). Hence, in principle, this state offers 
a convenient tool to follow the smooth transition from the particle to the wave 
nature of light. We have shown that expectation values computed in the case of 
an initial PACS at various instants during its temporal evolution in a Kerr-like 
medium offer a means to assess the effect of imperfect coherence of the initial state 
upon its subsequent dynamics p8|. 



A detailed analysis of the squeezing and higher-order squeezing properties of 
the wave packet in the neighborhood of a two-sub-packet fractional revival provides 
another method to ascertain whether the initial state is coherent or not |29|]. A 
quantifier to assess the extent of non-classicality of a state can be identified f30| 



in terms of the Wigner function corresponding to that state. We have used this 
quantifier to also identify the manner in which non-classical effects in the dynamics 
of the wave packet are affected, when the initial state is itself a non-classical state 
such as a PACS. The Wigner function is a quasi-probability distribution function 
of (a pair of) phase space variables, and could become negative in some regions 
of the phase plane. It is well known that the Wigner function corresponding to 
a CS (which is considered to be a classical state) is positive in all regions of the 
phase plane. This remains true for the Wigner function of the squeezed vacuum 
state as well. It is therefore not unreasonable to regard the negativity of this 
function in some region of the phase plane as an indicator of the non-classical 
nature of the state concerned. Further, the Wigner function can be reconstructed 
from experimentally measured quantities: for instance, that of a one-photon Fock 
state has been reconstructed using a phase-randomized pulsed optical homodyne 



method 31 . The indicator pOl based on the Wigner function that we have used 



is merely one way of quantifying the extent of the non-classicality of a state, but 
it suffices for our purposes. The Sudarshan-Glauber P-function [P^j , for instance, 
would provide another measure of the non-classicality of a state. 

The dynamics of a single mode of the radiation field enables us to understand, 
to a considerable extent, the connections between the behavior of quantum expec- 
tation values and revival phenomena. When two-mode interactions are considered, 
another significant aspect of quantum physics appears — namely, entanglement. 
We have examined the conditions under which analogs of revivals and fractional re- 
vivals of the initial state occur in this case. Revivals of two initially non-entangled 
modes governed by a nonlinear Hamiltonian have been examined in the literature 
||33|| . However, in order to study the full spectrum of randomness exhibited by 
the expectation values of relevant observables as they evolve in time, we need a 
system in which revival phenomena can either occur or be suppressed, depending 
on the values of the parameters in the Hamiltonian considered. A good candidate 
Hamiltonian for this purpose is obtained by a straightforward extrapolation of the 
single-mode example considered earlier, to explicitly take into account the atomic 
modes in the Kerr-like medium through which a single-mode electromagnetic field 
propagates, and the interaction between the atoms and the field. The Hamiltonian 
for the total system is given by pi 



H = hu a} a + fujotfb + h^tf 2 b 2 + kg (a)b + 6 f a). (1.9) 

As before, the ladder operators a and a) pertain to the field, while b and b^ are the 
ladder operators for the nonlinear oscillator modeling an atom of the medium. 7 is 
the anharmonicity parameter, and g is the strength of the interaction between the 
field and the medium. It is easily verified that the operator (a^a + b^b) commutes 
with H. We have analyzed this model for initial states that are direct product 
states, in which the atom is in its ground state, while the field is, respectively, in a 
Fock state, a coherent state, and a photon-added coherent state. Detailed results 
have been obtained on the precise manner in which the moments of appropriate 
quadrature variables display the counterparts of the revival and fractional revival 
phenomena in this case. These phenomena arise for weak nonlinearity, i.e., for 



S 



sufficiently small (<C 1) values of the ratio j/g. 

We have shown that these revival phenomena also manifest themselves in the 
entropy of entanglement of the system. The latter may be taken to be either 
the sub-system von Neumann entropy or the sub-system linear entropy, where by 
"sub-system" we mean either the field or the medium. These entropies are defined 
in terms of the corresponding reduced density operators. Striking differences arise 
in the behaviour of the entropies for different initial states of the field. These have 
been described in detail in the thesis. 

The situation that prevails for sufficiently large values of j/g is especially inter- 



esting. We have carried out a time-series and power spectrum analysis |35|, [36], |37 
of the time series generated by the values of the mean photon number (af a) com- 
puted over long intervals of time, and assessed the effect of the departure from 
coherence of the initial field mode on the ergodicity properties of the mean photon 



number The time series analysis shows that a wide range of behavior can 

manifest itself in the dynamics of the sub-system represented by the field mode: 
a progression occurs, from mere ergodicity for a coherent initial state and weak 
nonlinearity, all the way to "chaotic" behavior (as characterized by a positive 
Lyapunov exponent) for a sufficiently large departure from coherence of the initial 
state (i. e., a PACS with a relatively large number of added photons) and a suffi- 
ciently large nonlinearity (as given by the ratio 7/(7). This is corroborated by the 
fact that similar behavior is exhibited by (b^b), in such a manner that (6 t 6) + (a t a> 
remains constant, as required. The implications of "chaotic" dynamics, and what 
it really means in this context, are of special interest. A discussion of these aspects 
is given in the relevant chapter of the thesis. 

A summary of the contents of the contents of the rest of this thesis is as follows: 

Chapter 2 is devoted to a discussion of the manner in which distinctive sig- 
natures of revivals and fractional revivals are displayed in the time evolution of 
appropriate observables, which help in selectively identifying different fractional 
revivals |I8|. The model Hamiltonian used is that of a single mode propagating 



in a Kerr-like medium, Eq. (jl.Sf) , while the wave packet considered corresponds 
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to a coherent initial state. 

Chapter 3 describes how revivals and fractional revivals of a wave packet prop- 
agating in the nonlinear medium provide signatures of the degree of coherence of 
the initial wave packet. The investigation focuses on revival phenomena in the 
case of an initial photon-added coherent state. A comparison vis-a-vis the cor- 
responding results obtained for an initial coherent state enables the elucidation 
of the role played by the coherence of the initial state in determining the revival 



properties of the state |28 



Chapter 4 is concerned with the squeezing and higher-order squeezing properties 
of a PACS as it propagates in the nonlinear medium. Once again, the precise 
manner in which departure from coherence of the initial state affects its squeezing 



properties at the instants of revivals and fractional revivals is brought out f29 



Chapter 5 is based on the Hamiltonian in Eq. (l.S), that describes the entangle- 
ment dynamics of a single-mode electromagnetic field interacting with a nonlinear 
medium, an atom of the latter being modelled as an anharmonic oscillator. Re- 
vival phenomena are now examined in terms of the expectation values of certain 
observables, as well as the sub-system von Neumann entropy and the sub-system 
linear entropy. The initial states considered are direct product states, in which the 
atom in its ground state, while the field is taken to be in either a number state, 
or a CS, or a PACS. Once again, this leads to an understanding of the role played 
by the initial state. 

Chapter 6 deals with the case in which revivals of the initial state are suppressed. 
The dynamical behaviour of the mean photon number is investigated |38[] by means 
of a time-series analysis, in order to identify the effects of different initial states on 
the entanglement and ergodicity properties of the system. This behavior is seen to 
range from quasi-periodicity to exponential instability, depending on the specific 
initial state and the values of certain parameters occurring in the Hamiltonian. 

Chapter 7 concludes the thesis with some brief remarks placing this work in a 
broader perspective, and a list of open problems and avenues for further research. 
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CHAPTER 2 



Manifestations of wave packet revivals in the 
moments of observables - I 



2.1 Introduction 

Revivals and fractional revivals of a wave packet propagating in a nonlinear 
medium have been extensively investigated in the literature ||. Since revival 
phenomena arise from quantum interference, it is of interest to examine precisely 
how their occurrence is captured in the expectation values of appropriate physical 
observables. Such investigations have led to a better understanding of different 
aspects of revival phenomena in a wide variety of physical situations, ranging from 
wave packet propagation in the infinite square well |39| to the long-time behav- 
ior of a quantum bouncer In this Chapter, we discuss the manner in which 



distinctive signatures of revivals and fractional revivals are displayed in the time 
evolution of certain observables, thereby facilitating the selective identification of 
different fractional revivals ]nj] . It is necessary to bear in mind the salient features 
of the specific quantum interference between the basis states comprising the wave 
packet when we analyze the dynamics of expectation values. For this purpose, 
we outline in the next Section the conditions under which revivals and fractional 
revivals of a wave packet occur. 



2.2 Revivals and fractional revivals in wave packet 
dynamics 

The broad features of wave packet dynamics are quite generic, regardless of the 
details of the physical system concerned, the initial wave packet considered, and 



the specific nonlinear Hamiltonian governing the time evolution. We recall from 
Chapter 1 that, given an initial state |^(0)), a revival is the return of the overlap 
function 

c(t) = \(m\m)\ 2 (2.i) 

to its initial value of unity at specific instants of time. Revivals and fractional 
revivals are essentially controlled by the parameters occurring, respectively, in the 
first and second-order terms in the Taylor expansion of the energy spectrum E n 
about the energy E m corresponding to the peak of the wave packet. Hence, it 
suffices to consider wave packet evolution governed by an effective Hamiltonian 
whose energy eigenvalues are at most quadratic functions of n. Further, one of 
our objectives is to examine the role played by the initial state in the subsequent 
dynamics. Therefore, for purposes of comparison, an appropriate reference state 
must be selected. A good starting point is the investigation of the revival phenom- 
ena displayed by an initial coherent state (CS) propagating in a Kerr-like medium, 
so as to facilitate subsequent comparison with the dynamics of initial states which 
depart in a quantifiable manner from coherence. We explicitly demonstrate below 



the role played by quantum interference in this context |L3 



The effective Hamiltonian for the propagation of a single-mode electromagnetic 
field in a Kerr medium is given by JIT], ^ 



iJ = Xa t2 a 2 = X N(N-l) (2.2) 

in the usual notation, where a (respectively, a) ) is the photon annihilation (respec- 
tively, creation) operator, N = a)a is the photon number operator, and x (> 0) is 
essentially the third-order nonlinear susceptibility of the medium. Here, and in the 
rest of this thesis, h has been set equal to unity. We shall restore this factor when 
it becomes essential to do so (e. g, in Sec. 6.1 of Chapter 6). This Hamiltonian is 
also relevant in a very different physical context, namely, in describing the dynam- 
ics of a Bose-Einstein condensate in a potential well. In that case, a and a) are 
boson annihilation and creation operators, and \ characterizes the energy needed 
to overcome the inter-atomic repulsion in adding an atom to the population of the 
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potential well ||43|| . Many of the results that we derive may therefore be expected 



to be applicable in the context of atom optics and Bose-Einstein condensate, as 
well. 

The initial state of the field \ip(0)) is taken to be the CS \a), where a is any 
complex number. It satisfies 

a \a) — a \a) . (2.3) 

The Fock state representation of the CS, which we will use extensively, is given 
by 

| a) = e -M 2 /2£ « | n> . (2 . 4) 



n=0 v'n! 



The unitary time evolution operator U(t) corresponding to the Hamiltonian of Eq. 



is 

U(t) =exp[-ix*N(N -1)]. (2.5) 
Hence the state of the field at time t is given by 

|V(*))=exp[-ix*N(N-l)]|a). (2.6) 

The time evolution operator at instants t = 7r/(kx) (where k is an integer), namely, 

?7(7r/A;x)=ex P [--N(N-l)], (2.7) 

displays interesting periodicity properties. These follow from the fact that the 

eigenvalues of N are integers. We will use N to denote these eigenvalues. For odd 
integer values of k, 

exp [ — — (iV + k)(N + k - 1)] = exp [^-iV(iV - 1)]. (2.8) 

Similarly, for even integer values of fc, 

1 TT ITT 

exp [-- -(N + k) 2 } = exp [~jN 2 ]. (2.9) 
As a consequence, U(n/kx) can be expanded in a Fourier series with exp (—2nij/k) 
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as the basis functions, in the form 



exp [-jN(N - 1)] = |^exp[-^iV] (2.10) 

3=0 



for odd values of k, and 



exp[-fiV 2 ] = t>exp[-^iV] (2.11) 



j=0 



for even values of k, respectively, where the coefficients fj and gj are known. Using 
the above equations and the property 

e l ^ a \a) = \ae ix ), (2.12) 

we find 

f Z3S/i|«e-*^/*>, fc odd; 

|-0(7r/fc X )> = < (2-13) 
I Ejifij lae^e- 2 "*^*^ fc even. 

It is evident that at time T rev = 7r/%, corresponding to A; = 1, the initial 
state revives for the first time. Periodic revivals occur at integer multiples of 
T rev . Further, between t = and t = T rev and between two successive revivals, 
fractional revivals occur at instants nT rev + T rev /k. At these instants, the initial 
wave packet evolves to a state that can be described as a finite superposition of 
'rotated' coherent states with definite amplitudes. For instance, the state at time 
ir/(2x) (corresponding to k = 2) is a superposition of the two coherent states 
\ia) and |— ia). In general, \ip(ir/kx)) is a superposition of k coherent states. 
The corresponding wave packet in position space is a superposition of k spatially 
distributed Gaussian wave packets. The periodicity property of U further implies 
that the wave function at times t = irj/ (kx), 1 < J < A: — 1 for a given k, is also 
a superposition of k wave packets [011 . 
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2.3 Signatures of revivals and fractional revivals 
in quantum expectation values 



We now proceed to examine the manner in which fractional revivals are manifested 
in distinctive ways in the expectation values of the physical observables pertaining 
to the system at hand. As the system enjoys revivals with a period T rev , all such 
expectation values are periodic functions of t with this fundamental period. A 
family of relevant observables is provided (in the context of nonlinear media) by the 
field quadrature (a e tip + a^ e~ tlf ) for various values of the phase ip. In the context of 
Bose-Einstein condensation, the expectation value (a(t)) of the atom annihilation 
operator represents the condensate wave function. Its real and imaginary parts 
(which correspond to the cases (p — and ip = — |tt, respectively) are analogous 
to classical phase space variables. Let us define the operators 

(a + ah , (a — ah . 
x = - ■=-!- and p = - =-A 2.14 

Clearly, their expectation values alone do not suffice to reproduce the full informa- 
tion contained in the wave function itself. In principle, an infinite set of moments, 
comprising the expectation values of all powers of x and p and their combinations, 
is required for this purpose. In this sense, the quantum system is equivalent to an 
infinite-dimensional classical system in which the role of the dynamical variables 
is played by the set of expectation values. However, we emphasize that even the 
first few moments can be shown to yield considerable information on the behavior 
of the system. 

Recalling that a \ a) — a \ a), we define the c- number function 

a{t) = {i/>(t)\a\il>{t)) = {a\e iHt ' h ae- iHt ' h \a), (2.15) 
so that a(0) = a. For the case at hand, this simplifies after some algebra to 
a(t) = ae ~H 2 (i~cos2 X t) ^ cog Q a j2 sin (2 x t)) -isin (H 2 sin(2xt))] . (2.16) 
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Thus a(t) is a periodic function of time with period ir/x, as expected. It is 
convenient to introduce the notation 



and 



a = a\ + 1012 = -1= (2.1 7 ) 

a/2 



l«| 2 = ^o+Po)- (2-18) 



Xo and po represent the locations of the centers of the initial Gaussian wave packets 
corresponding to the CS \a) in position space and momentum space, respectively. 
The expectation values of x and p at any time can then be obtained as explicit 
functions of t (as stated in Appendix B) in the form 

(x(t)) = e -"(i-«*2x9 [ Xq cos(z/sin2xt) + Po sin (v sin 2xt)] , (2.19) 



(p(t)) = e -"(i-co S 2xt) [_ Xq S i n (^ S in2x^) +Pq cos {y sin 2 X t)] . (2.20) 

Expressions for the higher moments can be deduced readily from the general 
result 

(a tr a r+s ) = a s v r e - v ^- co& 2sxt) exp [-ix(s(a - 1) + 2rs) t - iv sin 2s X t] , 

(2.21) 

where r and s are non-negative integers. This result is derived in Appendix A 
(see Eq. (| A . 1 5| ) ) . It is a special case of a more general relation derived there (Eq. 



( |A.14Q ). Using this result, the second moments of x and p are found to be given 
by 

2 (x 2 {t)) = 1 + xl + p 2 Q + e ~ u ^- cos 4 *' } [{xl - pi) cos (2 X t + v sin A X t) 

+ 2x p sin (2xt + v sin 4 X t)] , (2.22) 

2{p 2 (t)) = l + xl+p 2 -e- u{1 - cosixt) [{xI-pD cos(2xt + z/sin 4 X t) 

+ 2x p sin (2%t + v sin 4^t)] . (2.23) 
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For reasons already mentioned, the higher moments also carry much information 
of direct interest. We give below the expressions for the third and fourth moments 
as well, since we need them to obtain the skewness and kurtosis in all the cases of 
interest. The third moments can be written compactly in the form 

4 (a; 3 (t) ) = e-^- cos 6 ^ [(a; 3 - 3x p 2 ) cos (6 X t + v sin 6x0 
+ (3xqP - pi) sin (6xt + v sin 6%*)] 

+ 6u [ (x(t)) (1 + cos 2x0 + (p(0> sin 2 X t)] , (2.24) 



4 (p 3 (t)) = e -(i-cos6 X t) [^3 _ 3a , op ^ gin ^ + ^^6^) 
— (3xqPo — Po) cos (^X^ + v si n 6x0] 

+ 6u [(p(f)) (1 + cos2x0 - (x(t)) sin 2#]. (2.25) 
The fourth moments are given by 

8<a; 4 (0) = e-^ 08 8 ^ [(4+p 4 -6x^) cos (I2 x t + Z / S in 8x0 

+ 4(xqPo - xoPo) sin (12x£ + v sin 8x0] 

+ 8z/e-" (1 - cos 4xt) [(x 2 - pi) cos (6 X t + v sin 4 X 

+ 2x £>o si 11 (6x^ + v si n 4x0] 

+ 24<x 2 (0>+3(^+^)2_ 6) ( 2 ,2 6 ) 



8<p 4 (0> = e-^- cos8 ^[(x 4 +^-6^)cos(12xt + z/sin8x0 

+ 4(xqPo - ^oPo) sin ( 12 x^ + v sin 8x0] 

- 8ue- u{1 - cos Axt) [(x 2 - pg) cos (Q X t + i/ sin 4x0 
+ 2x p sin (6x* + v sin 4x0] 

- 2A(x 2 (t))+3(xl+pl) 2 + 24(x 2 +pl) + 18. (2.27) 

With these expressions, the variances of x and p as functions of t, as also the 
skewness and kurtosis in each case, can be obtained. The uncertainty product 
(i. e., the product of the respective standard deviations) AxAp, which initially 
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has the minimum value ^ , is of special interest . We do not write down the lengthy 
expressions for these quantities here, but we shall comment upon their time vari- 
ation in the sequel, when discussing the numerical results for specific values of the 
parameters. 



2.4 Moments of x and p: A classical interpreta- 
tion 

Before discussing our results in detail, we digress briefly to compare the explicit 
solutions found for the expectation values of x and p in Eqs. ( |2.19| ) and ( |2.20| ) 



with the solutions that would have been obtained for x(t) and p(t) had the system 
been a classical one, governed by the classical counterpart of the normal-ordered 
Hamiltonian (restoring the factor h) H = hxa) 2 a 2 , namely, 

H cl = ^(x 2 +p 2 ) 2 . (2.28) 

Although the equations of motion corresponding to H c \ are nonlinear, it is evident 
that x 2 + p 2 is a constant of the motion, so that the phase trajectories are circles. 
However, the frequency of motion is dependent on the initial conditions (i.e., the 
amplitude of the motion), being equal to v — \{x\ + pi). This is, of course, a 
well-known feature of nonlinear oscillators. But we note that the actual solutions 



for (x(t)) and (p(t)) in Eqs. ( |2.19| ) and fl2.20|) are more complicated than the 



classical ones for x(t) and p(t) under H c \. This is a consequence of the quantum 
mechanical nature of the system, over and above the nonlinearity of H. However, 
the expressions for (x(t)) and (p(t)) can be given the following interesting inter- 
pretation in classical terms. Define the (non-canonical) pair of classical dynamical 
variables 

X = (x) e u(1 - cos2xt) , P=(p) e v{1 - cos2xt \ (2.29) 

and the re-parametrized time t = sin (2%t). The initial values of these variables are 
again x and po, respectively. Equations (|2.19| ) and (|2.20D can then be re- written 
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in the suggestive form 

X = x cos vr + po sin vr , P = — xq sin vr + po cos vr. (2.30) 

But these are the solutions to the system of equations 

dX/dr = vP , dPjdr = -uX, (2.31) 

describing a nonlinear oscillator of frequency 

v = \(xl+pl) ) (2.32) 

in terms of the transformed variables (X, P) and the re-parametrized time r. At 
the level of the first moments, therefore, the system is effectively a nonlinear 
oscillator after a suitable transformation of the relevant variables, together with a 
re-parametrization of the time. 



2.5 Signatures of fractional revivals 

We now discuss our results, using representative numerical values for the vari- 
ous parameters in the problem. Turning to the time dependence of the various 
moments of x and p, there are two striking features that underlie the essential 
point we wish to make. First, the higher the order of the moment (or cumulant), 
the more rapid is its variation, since the leading (highest) frequency in the k th 
moment is 2k\- Second, the time dependence is strongly controlled by the factor 
exp [—1/ (1 — cos 2kxt)], k — 1,2,..., that modulates the oscillatory terms. While 
this permits substantial time variation for sufficiently small values of u, it acts 
as a strong damping factor for large values of z/, except when cos (2mxt) is near 
unity. As one might expect, this happens precisely at revivals (when t = nn/x, 
an integer multiple of T rev ). But it also happens — in the k th moment and not in 
the lower moments — at the fractional revival times t = (n + j /k)T rcv . Thus, by 
setting v at a suitably large value, we can ensure that the moments are essentially 
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fa) (b) (c) 

Figure 2.1: (x) as a function of time 



static, bursting into rapid variation at specific instants of time before reverting to 
quiescence. 

These points are illustrated in the figures. Owing to an obvious symmetry of 
H, the moments of x and p behave in an essentially similar manner, especially 
if we start with the symmetric initial condition xq = po- Without significant 
loss of generality, we restrict ourselves to this case in what follows. We have set 
X = 5 in the numerical results to be presented, but this is irrelevant, as all the 
plots correspond to t measured in units of T rev = tt/x- We find that for very 
small values (<^C 1) of xq and po (i. e., of u), the nonlinearity of H does not play 
a significant role, and the behavior of the system is much like that of a simple 
oscillator. Interesting behavior occurs for larger values of v. We therefore present 
results for three typical values of the parameters representing the initial conditions, 
namely: 

(a) xq = po = 1 v = 1; 

(b) x = po = a/10 =^ v = 10; 

(c) x = Po = 10 =^ v = 100. 

These correspond, respectively, to small, intermediate, and large values of v. In 
all the "phase plots", the point representing the state at t = is labeled A. 

Figures |2.1| (a)-(c) show the variation of (x(t)) as a function of t for small, 
medium and large values of v. (As already mentioned, (pit)) displays similar 
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(a) (b) (c) 

Figure 2.2: "Phase plot" of (p) versus (x) 

behavior.) For sufficiently large values of u, it is evident that, except for times 
close to integer multiples of T rev , (x(t)) and (p(t)) essentially remain static at 
the value zero. Figures ^2|(a)-(c) depict the corresponding "phase plots" in the 
((x) , (p)) plane. In Fig. |2T^(c) , the representative point remains at the origin most 
of the time, except at times close to successive revivals, when it rapidly traverses 
the rest of the curve before returning to the origin. 

While sudden changes from nearly static values of (x(t)) and (p(t)) are thus 
signatures of revivals, the occurrence of fractional revivals is not captured in these 
mean values. The fractional revival occurring mid-way between successive revivals 
(e. g., at t = tc/2x in the interval between t = and t = T rev ), when the initial 
wave packet reconstitutes itself into two separate wave packets of a similar kind, 
leaves its signature upon the second moments. Figures |2.3| (a)-(c) show the vari- 
ation with time of the uncertainty product Ax Ap. In each case, this product 
returns at every revival to its initial, minimum, value (= |), rising to higher val- 
ues in between revivals. Once again, for sufficiently large values of z/, the product 
remains essentially static at the approximate value (~ + v) for most of the time, 
but undergoes extremely rapid variation near revivals, and also near the fractional 
revivals occurring mid-way between revivals. During the latter, the uncertainty 
product drops to smaller values, but does not reach the minimum value |. 

There is a very striking difference in the behavior of the standard deviations 
near revivals as opposed to their behavior near the foregoing fractional revivals. 
This is brought out in Figs. |2.4j (a)-(c), which is a "phase plot" of Ap versus 



21 




Oh 
< 

< 



m 



v=10 



< 60 
X 

< 

-40 



V=1QQ 



t/T 



t/T r 



(a) (b) (c) 

Figure 2.3: Variation of the uncertainty product with time 
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Figure 2.4: "Phase plot" of Ap versus Ax for v — 1, 10 and 100, respectively. 

Ac. For very small i/, as in Fig. |2.4j (a), Z\x and Z\p vary quite gently around a 
simple closed curve. When v is somewhat larger, as in |2.4| (b) which corresponds 
to v = 10, the plot begins to show interesting structure. For much larger values of 
v as in |2.4j (c), the initial point A quickly moves out on the zig-zag path about the 
radial Ap = Ax line to the steady value represented by the point B, and returns to 
A at every revival along the complementary zig-zag path. Close to the fractional 
revival at t = (n + |)T rev , however, the representative point moves back and forth 
along the azimuthal path BCDB rather than the zig-zag path: clearly, a kind of 
"squeezing" occurs, as one of the variances reaches a small value while the other 
becomes large, and vice versa. (Of course the state of the system is far from a 
minimum uncertainty state throughout, except at the instants nT rev .) 

The fractional revivals occurring at t = (n + |)T rev and t = (n + |)T rev , at 
which the initial wave packet is reconstituted into a superposition of three separate 
wave packets, are detectable in the third moments of x and p. To make this 
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Figure 2.5: Square of the skewness in x as a function of time 




Figure 2.6: The "phase plot" of (3^ versus (3\ 

unambiguous, we may consider the third moments about the mean values — or, 
in standard statistical notation, the square of the skewness, defined as 



0\ 



(x) 



((x 2 ) - (x) 2 ) 3 ~ (Ac) 



(2.33) 



and similarly for (3[ p \ Figures p.5| (a)-(c) show the variation of (3^ with t. It is 
evident that, for sufficiently large values of v, (3^ remains nearly zero most of the 
time, except for bursts of rapid variation close to revivals and fractional revivals. 
Both fii and (3^ actually vanish at t = nT mv (since the Gaussian wave packet 
corresponding to the CS \a) has no skewness in either x or p), but they remain 
non-zero at t — (n+ |j) T rev , j = 1,2. More detailed information is obtained from 
a "phase plot" of /3^ versus /3[ x \ depicted in Figs. (a)-(c). 



Finally, we consider fractional revivals corresponding to m = 4, when four 
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superposed wave packets appear. These are detectable in the behavior of the 
fourth moments of x and p. Equivalently, we may use the excess of kurtosis 
ifli — 3), where the kurtosis of x is defined as 



((x 2 ) - (x) 2 ) 2 ~ (Ac) 



^4 ' 



(2.34) 



with a similar definition for The excess of kurtosis is the measure of the depar- 
ture of a distribution from a Gaussian. Figures |3(a)-(c) depict how ((3 { 2 x) - 3) 



varies with time. For sufficiently large v, both (/^ — 3) and ((3^' — 3) re 
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Figure 2.7: The excess of kurtosis of function of time 
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Figure 2.8: "Phase plot" of {$> - 3) versus (f3f> - 3) 

main essentially static near the value — | most of the time. They vary rapidly 
near revivals, vanishing at t = nT rev because the wave packet is a Gaussian 
both in position space and in momentum space at these instants of time. As 
is clear from Fig. |2]7](c), they also vary rapidly near the fractional revivals at 
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t = (n + j/4)T rev (where j = 1,2,3), oscillating about the "steady value" — |. 
Once again, a "phase plot" of {(3^ — 3) versus (/Jj — 3) helps identify features 
that distinguish between the three fractional revivals concerned. These are shown 
in Figs. [2]^ (a)-(c). 

2.6 Concluding remarks 

We have shown that distinctive signatures of the different fractional revivals of a 
suitably prepared initial wave packet are displayed in the mean values and higher 
moments of appropriate observables. The complicated quantum interference ef- 
fects that lead to fractional revivals can thus be captured in the dynamics of these 
expectation values, which may be regarded as effective dynamical variables in a 
classical phase space. While direct experimental measurement of the expectation 
values considered, particularly the higher moments of the quadrature variables, is 
certainly not easy, it is not, a priori, impossible. Moreover, the essential point 
of our analysis is to bring out explicitly an interesting aspect of quantum inter- 
ference: the link between the occurrence of a fc-sub-packet fractional revival and 
the behavior of the fcth moment of appropriate observables, emphasizing (as in 
the rest of this thesis) the role of quantum expectation values in the investigation 
of non-classical effects. Again, it is evident that the specific observables signaling 
the occurrence of revivals and fractional revivals will vary from system to system, 
and it may not be easy, in practice, to identify the correct observables in all cases. 
However, the special nature of the interference between basis states at certain 
definite instants leads to revival phenomena, and this feature is closely linked to 
sudden changes in these expectation values at those instants. This feature is not 
necessarily traceable in all cases to the occurrence of an overall exponential factor 
in the expressions for the expectation values, although our analysis seems to in- 
dicate that a factor involving an exponential of periodic functions certainly plays 
a crucial role. Wave packet propagation in a Kerr-like medium provides us with 
a clean framework for illustrating these aspects, as it involves just a single pair 
of ladder operators, and the spectrum is not degenerate, in general. In the next 
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Chapter we therefore extend our investigations to initial photon-added coherent 
states propagating in a Kerr-like medium, and draw attention to the differences 
that arise due to the departure from perfect coherence of the initial state. 
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CHAPTER 3 



Manifestations of wave packet revivals in the 
moments of observables - II 



3.1 Introduction 

In this chapter, we study the manner in which signatures of the revival phenom- 
ena which are manifested in the expectation values of observables depend on the 
extent of coherence enjoyed by the initial wave packet. For this purpose, we 
need to compare the details of the revival phenomena exhibited by an appropriate 
non-Gaussian initial state with that of an initial CS propagating in the nonlinear 
medium. As mentioned in Chapter 1, a suitable candidate for the initial state is 
the photon-added coherent state (PACS): it possesses the useful property of a pre- 
cisely quantifiable, as well as tunable, degree of departure from perfect coherence. 
This property follows from the fact that a PACS can be obtained in principle by 
repeated addition of photons to a CS. The normalized m-photon- added coherent 
state \a,m) is defined as HI] 



(qt)"|q) (at)™ !<») 

a, m) = — . = — . . (3.1) 

A/(a|a m at m \a) y/m\L m (-v) 

where m is a positive integer, v = \a\ 2 as before, and L m (—u) is the Laguerre 
polynomial of order m. Setting m = 0, we retrieve the CS |a). The mean photon 
number in the PACS \a, m) is given by |41 



(m+DW-,) ^ (3 . 2) 



When v = 0, (N) m reduces to m, as required. When v m, we can show that 
(N) m = v + 2m + 0{y~ x ). The plot of ((N) — v) versus v in Fig. |3.1j shows how 
the ordinate saturates to the value 2m as v increases. The extent of departure of 




a PACS from perfect coherence clearly becomes more pronounced with increasing 
m. A PACS also possesses other interesting properties which may be kept in mind 
when we investigate its dynamics. For instance, in contrast to a CS, it exhibits 



phase-squeezing and sub-Poissonian statistics ]T3[ . The latter implies that the 

standard deviation of the photon number operator N behaves like (N) 2 rather 
i 

than (N) 2 , the exponent (3 being a calculable decreasing function of m. 

The PACS is also a nonlinear coherent state in the sense that it is an eigenstate 
of a nonlinear annihilation operator [ill, namely, 



m \ 

1 — a\a,m) = a \a, m) . (3.3) 

1 + a* a J 

The state \a, m) can also be viewed in yet another way. Instead of the Fock basis 
{|n}}, we may consider the unitarily transformed basis {\n, a)} formed by the 

generalized coherent states 

\n } a) = e aaWa \n) . (3.4) 

Equivalently, for a given n, the generalized coherent state |n, a) is simply the state 
(g^ — a*) n \ot), normalized to unity. It can now be shown that, for a given m, the 
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Figure 3.2: Probability of single-photon content versus v in a CS (m = 0) and in 
a 1-photon- added CS (m = 1). 



state \a, m) is a finite superposition of the form J2n=o Cn l n > a )- 

We further note that the single-photon-added coherent state recently generated 



in experiments |27| has a higher probability pi of single-photon content than an 
ideal CS for sufficiently small values of v, as may be seen from Fig. |3.2j . (It is 
straightforward to show that this happens for v < (vo — l)/2 ~ 0.618.) In this 
sense, a PACS would seem to be preferable to a CS for single-photon experiments. 



With this brief introduction to the properties of a PACS, we proceed to examine 
the details of its dynamics as it propagates in a nonlinear medium. The latter is 
once again modeled by the Hamiltonian of Eq. (|2.2|) . 



3.2 Dynamics of photon- added coherent states 
propagating in a nonlinear medium 

As in the case of an initial CS, we examine the moments and variances of the 
hermitian operators x and p, defined in Eq. ( [2.14| ), as the single-mode field evolves 



from the initial PACS of Eq. (3J.) under the Kerr-like Hamiltonian. We shall use 
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the convenient notation 



(x(t)) m = (a,m\ e m l n xe- im l % \a , m) , (3.5) 

with an analogous definition for (p(t)) m . As H is diagonal in the number operator 
N, it follows that the mean photon number (N) m , as well as its higher moments, 
and hence the sub-Poissonian statistics of the photon number, remain unaltered 
in time. 

We recall from Chapter 2 that, at the level of the first moments, the dynamical 
equations in the case of an initial CS are those of a classical nonlinear oscillator, 
with a certain re-parametrization of time. In contrast, the time dependences of 
{x(t)) m and (p(t)) m (i. e., averages obtained in the time-evolved m-photon-added 
coherent state with m 7^ 0) differ in striking ways from the foregoing, even for 
small values of m. The revival time T rev remains equal to tt/x, of course, but 
in the intervals between revivals the time evolution is considerably more involved 



than the expressions in Eqs. Q2.19Q and fl2.20[) for the case m — 0. This implies 
that even a small departure from coherence in the initial state and from Poissonian 
number statistics leads to a very different time evolution of the system and the 
phase squeezing it exhibits. We have calculated analytically the exact expressions 
for (x(t)) and (p(t)) in Appendix B, and these are given in Eqs. ( B.3|) and 
( |B.4|) . The initial values of these expectation values in the PACS \a,m) are 

<*(0)>m = J^TV, X « ' (P(°)>m = TT^ ' ( 3 ' 6 ) 

where L l m {— v) = d L m+1 (— v)jdv is an associated Laguerre polynomial, x = 
y/2Kea, and po = \/2lma. Analogous to the variables X and P in the case of 
the CS given by Eq. (|2~29|) of Chapter 2, let us define 



X m {t) = (x(t)) m exp [v (1 - cos2x*)], 

Pm(t) = ip{t))m eX P I" (! " COs2 X t)]. 



(3.7) 



(Thus Xq and Pq are just the quantities X and P of the preceding chapter.) The 
solutions for these quantities may then be written in the following compact and 
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suggestive form (derived in Appendix B, Eqs. (|B.11| ) and Eqs. (|B.12|) ): 



X m (t) = x Re z m (t) + p Im z m (t) 
P m (t) = Po Re z m (t) - x Im z m (t) , 



(3.8) 



where 



Li 




exp [i (2mxt + v sin 2%£)] . 



(3.9) 



A number of differences between these results and those for the case m = are 
noteworthy. First, \z m \ varies with t, in contrast to |2o(t)| = 1. (Note also that 
z m (0) = L m (—v) / L m {—v) 7^ 1.) The time dependence of X m and P m involves the 
sines and cosines of the set of arguments (2xlt + v sin 2%t), where I — m, . . . , 2m. 
Thus, not only are higher harmonics present, but the arguments also involve secu- 
lar (linear) terms in t added to the original [y sin 2\t.). This important difference 
precludes the possibility of subsuming the time dependence into that of an effec- 
tive nonlinear oscillator by means of a re-parametrization of the time, in contrast 
to the case m — 0. 

We present the rest of our results with the help of figures based on numerical 
computation. As in Chapter 2, we have set \ = 5 for definiteness. (Our results 
on wave packet dynamics will hold good for any other value of \ as well. The 
numerical value of x merely sets a time scale.) Again, as before, we restrict 
ourselves to the case xq = po (there is no significant loss of generality as a result 
of this symmetric choice of parameters). As before, the presence of the overall 
factor exp [—1/(1 — cos2xt)] in the expressions for (x(t)) and (p(t)) m implies 
that, for sufficiently large values of the parameter u, the expectation values remain 
essentially static around the value zero, and burst into rapid variation only in the 
neighborhood of revivals. Smaller values of v enable us to resolve the details of 
the time variation more clearly. 

Figures |0| (a) and (b) are, respectively, plots of the expectation values {x(t)) 1 
and (x(t)) 10 versus t/T rev for parameter values x = p = 1 (i.e., for v = 1). The 
revivals at integer values of t/T rev are manifest. With increasing m (or a decreasing 
degree of coherence in the initial state), the relatively smooth behavior of (x(t)) 
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t/T. 



Figure 3.3: (x(t)) as a function of t/T Tev for a PACS (a) \a, 1) and (b) \a, 10), with 
v = \. 



gives way to increasingly rapid oscillatory behavior in the vicinity of revivals. This 
is to be compared with Fig. |2.1| (a), the corresponding plot for an initial CS. The 
range over which the expectation value varies also increases for larger values of m. 

Essentially the same sort of behavior is shown by (p{t)) m . However, a "phase 
plot" of (p(t)) m versus (x(t)) m in Figs. |3]4] (a) and (b) reveals complementary 
aspects of such oscillatory behavior with increasing m, showing how the oscillations 
in the two quantities go in and out of phase with each other. The entire closed 
curve in each case is traversed in a time period T rev . 

As in the case of an initial CS, here too the initial state undergoes fractional 
revivals in the interval between any two successive revivals at the instants T rev /k 
where k = 2,3,... , and also, for any given k, at the instants jT rev /k where 
j = 1, 2, . . . , (k — 1). These fractional revivals at the instants j T rev /k also show 
up in the rapid pulsed variation of the k th moments of x and p, and not in the 
lower moments. However, there is a crucial difference if we use \a, m) as the initial 
state: then, even for relatively small values of m, the signatures of fractional 
revivals appear for values of v that are not large, in contrast to what happens 
when the initial state is the coherent state |a) . (Recall that (N) m is determined 
by v according to Eq. (|3.2|).) 



For illustrative purposes we investigate the specific case of the fractional revival 
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10 

Figure 3.4: "Phase plot" of (p) vs (x) for an initial PACS (a) \a, 1) and (b) \a, 10), 
with v = \. 




Figure 3.5: Phase distribution P(4>) versus <p in units of tt, for v = 0.3 and m — 0, 1 
and 10. 



at t = 7|T rev . This corresponds to the appearance of two spatially separated similar 
wave packets, i. e., a single qubit in the language of logic gate operations. 

At this instant the overlap between the two sub-packets in the wave function 
corresponding to the state |^(T rev /2)) decreases significantly with an increase in 
m. This feature emerges clearly in a plot of the phase distribution P(<f>) versus 
the angular variable <fi (shown in Fig. [T5] for v = 0.3 and m — 0, 1 and 10). This 
phase distribution is defined as |45l 

1 oo 

P(<P) = -\(<PmTr ev /2))\ 2 , |0>=5>^|n). (3.10) 

n=0 

In this sense, a PACS is better suited than a CS for performing single qubit oper- 
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Figure 3.6: Variation of the uncertainty product with t/T rev for an initial CS 
(m = 0) and initial PACS \ce, 1) and \a, 10) with v — 1. 

ations treating the two sub-packets as qubits |0) and |1), respectively. Similarly, 
at t = |T rev , the overlap between the four sub-packets is significantly smaller in 
the case of an initial PACS than in the case of an initial CS, and this feature be- 
comes more pronounced with increasing m. The sub-packets therefore behave as 
genuine orthogonal basis states in the former case even for relatively small values 
of v. Hence the quantum logic approach to wave packet control for the purpose 
of implementing two-qubit operations, outlined in Ref. ||25|| , is more feasible with 
an initial state that is a PACS rather than a CS. 

Plots of the product Ax Ap of the standard deviations of x and p versus t/T Tev 
over a full cycle are shown in Fig. [T6] for initial states given, respectively, by the 
coherent state |a) (dotted curve), the photon-added state \a, 1) (dashed curve), 
and the 10-photon-added state |a,10) (bold curve), for x — 5, x = p Q = 1, so 
that v — 1. It is seen that hardly any trace of the fractional revival is evident 
in the case of an initially coherent state, in marked contrast to the case of the 
photon-added states, in which the fractional revival is signaled by oscillations 
whose frequency and amplitude increase quite rapidly with increasing m. This 
effect gets masked for larger values of the parameter z/, when these oscillations are 
relatively insensitive to the value of m. Another striking feature that provides a 
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clear distinction between the revivals at t — nT rev and the fractional revivals at 
t — (n + 5) 2~rcv (where n is a positive integer) is illustrated in Fig. pT7| , which is 
a plot of Ap versus Ax . The dotted and full lines in the same plot correspond 
to m = and m = 5, respectively, making ready comparison possible. We recall 
from Chapter 2, Figs. |2.4| (a)-(c) that for an initial CS these standard deviations 
start from an initial value |. As t increases, they rapidly build up, oscillating 
about the radial Ap = Ax line with an initially increasing, and then decreasing, 
amplitude. A maximum value of Ax and Ap is attained, at which these quantities 
then remain nearly static, till the onset of the fractional revival at T rev . They 
then begin to oscillate rapidly once again, but this time in a tangential direction, 
swinging back and forth along an arc with an amplitude that initially increases 
and then decreases to zero: in other words, the individual standard deviations 
fluctuate rapidly in the vicinity of the fractional revival (while [(Ar) 2 + (Ap) 2 ] 1 / 2 
remains essentially unchanged in magnitude), in marked contrast to what happens 
at a revival. It is evident from Fig. [3.7| that all these features are very significantly 
enhanced and magnified for non-zero values of m, relative to what happens for 
the case m = 0. 

Similar signatures of fractional revivals for higher values of k can be picked out 
by using initial states \ot, m), even with relatively small values of v. The oscillations 
in the moments of observables become more pronounced with increasing m. For 
instance, signatures of the fractional revivals at t — \j T rev , j = 1, 2, 3 are clearly 
discernible in the behavior of the fourth moments of x. Figure |3.8| is a plot of {(Sp) 4 ) 
versus ((Sx) 4 ) where 5x = x — (x) , 8p = p — (p) , for v = 1. The dotted and full 
lines correspond to the initial states \a) and |a,2), respectively. 

3.3 Concluding remarks 

We have shown that distinctive signatures of wave packet revivals and fractional 
revivals are displayed in the temporal evolution of appropriate variables, enabling 
identification of departure from coherence of the initial state. These distinctions 
between an initial CS and an initial PACS are, however, clearly seen only for small 
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Figure 3.8: "Phase plots" of ((5p) 4 } and ((5x) 4 } for v = 1, for an initial CS and 
an initial PACS \a,2). 
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values of v. This is related to the fact that, for small values of v, the probability of 
the m-photon content of an m-photon-added coherent state is substantially larger 
than that of the CS. Another difference between these states emerges during their 
temporal evolution: the addition of even a single photon washes out the squeezing 
property displayed by an initial CS at specific instants. In the next Chapter, we 
examine the squeezing and higher-order squeezing properties of these states as 
they propagate in a Kerr-like medium. 
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CHAPTER 4 



Squeezing and higher-order squeezing in wave 

packet dynamics 



4.1 Introduction 

In this chapter, we investigate the squeezing and higher-order squeezing proper- 
ties of an initial PACS as it propagates in a Kerr-like medium. In particular, 
we examine the precise manner in which departure from coherence of the initial 
state affects its squeezing properties in the neighborhood of revivals and fractional 
revivals. 

For ready reference, we first summarize the conditions under which a generic 



state is said to display squeezing or higher-order squeezing | 46[] . Consider two 
arbitrary Hermitian operators A and B with a commutator 

[A,B]=iC. (4.1) 

The generalized uncertainty relation satisfied by these operators in any state is 

AAAB>^\(C)\. (4.2) 
The state is said to be squeezed in the observable A if 

{AAf < l -\(C)\. (4.3) 

A similar definition holds for B. 

In particular, for the quadrature variables x and p defined in Eq. ( f2.14j ), this 
implies that x (or p) is squeezed if Ax (or Ap) < 1/V2 in the state concerned. 



Several states satisfying this squeezing property have been investigated in the liter- 
ature || . Some examples are the squeezed vacuum state, squeezed number states, 
and photon-added coherent states. Proposals for the generation, and detection of 
squeezed states, as well as their potential applications, have been outlined in the 
literature [g7|, H • 



Of direct relevance to us is the squeezing property displayed by the m- photon- 
added coherent state \a,m). A general result for this state is as follows. Consider 
the hermitian operators 

06* + ate-** , „ ae i *-ate-*P , , 

satisfying the commutation relation [A v , B^] = i for any value of ip. For ip = 0, 
and B^ reduce to quadrature variables x and p. It can be shown that the variance 
of Ap in the state \a,m) is given by |41 



(^,) 2 = {2v(L* m (-v)L m (-v) - [Ll(-v)} 2 ) cos {2(6 + V )) - 2v[Ll n (-v)} 2 
- [L m (-z/)]V2(m+l)L m+1 (-z/)L m (-z/)| / /|2[L m (-z/)] 2 |, (4.5) 

where a = v l l 2 e ie as before. Setting m = in the above equation and simplifying, 
we retrieve the corresponding expression for a CS |a), namely, (AA^) 2 = |. As is 
well known, the CS is not a squeezed state. It is in fact a minimum uncertainty 
state in the observables x and p with Ax = Ap = l/y/2. Choosing the phases 
such that 9 + ip = 7r, a plot of (AA V ) 2 versus \a\ shows that almost 50% squeezing 
can be obtained over a wide range of values of |a|, for non-zero values of m. 

The concept of higher-order squeezing was first introduced by Hong and Man- 
del A state is said to be squeezed to order 2q (q = 1, 2, 3, . . .) in the operator 
A if ((5A) 2q ) in that state is less than the corresponding value obtained for a 
CS, where 5 A = A — (A). For the quadrature variable x, this reduces to the 
requirement 

<(^) 2 "><i(2g-l)!! (4.6) 
Another type of higher-order squeezing, called amplitude- squared squeezing, was 
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defined first by Hillery |50| , and generalized subsequently by Zhang et al. to q th - 
power amplitude-squeezing []5T|. (Amplitude-squared squeezing corresponds to the 
case q — 2.) One considers the two quadrature variables 



(a q + a) q ) , (a q -a) q ) . 

Zi = - — Z 7 ^r- L and z 2= —FT 1 (9 = 1,2,3,...)- 



V2 



iV2 



The generalized uncertainty principle now reads 



(4.7) 



{AZ x f{AZ 2 f>\\{[Z u Z 2 ])\\ 



The state is said to be q -power amplitude-squeezed in the variable Z\ if 



(AZ,) 2 < \ ( \Z X , Z 2 ] ) . 



(4.9) 



Amplitude squeezing in Z 2 is similarly defined. Let us write [a q , a) q ] = F q (N). 
This is a polynomial of order (q — 1) in the number operator N (see Appendix C, 
Eq. (lO! ), given by 



n=l ^ ' 



(4.10) 



Let 



Ar(*) 



l(f,(N)) 



(4.11) 



(The time-dependence has been indicated explicitly in D q (t) to remind us that the 
expectation values involved are those in the instantaneous state of the system.) It 
is easily seen that the state is g th -power amplitude-squeezed in Z\ if — 1 < D q < 0. 
We can rewrite Eq. (|4.11| ) in terms of a q and as 



D q (t) 



Re (a 2q ) - 2 (Re (a q ) ) 2 + {a) q a q ) 



(F q (U)) 



(4.12) 



We now examine the manner in which the extent of coherence of the initial 
state affects its subsequent squeezing properties. There is a considerable literature 
on this subject. For instance, dynamical squeezing of a PACS arising as a result 
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of the time-dependence of the frequency of the electromagnetic field oscillator in a 
cavity has been investigated ||52|| . The non-classical effects exhibited by an initial 
CS propagating in a Kerr-like medium |53| , its Hong-Mandel squeezing properties 



54| , and its q -power amplitude-squeezing properties [BHj have been discussed in 



the literature. We are interested specifically in the squeezing properties of wave 
packets propagating in a Kerr-like medium. In the next Section, we therefore 



summarize the results of |55j in this context, obtained for an initial CS. This will 
facilitate subsequent comparison between our work on the squeezing properties of 
an initial PACS [EjJ with those of an initial CS. 



4.2 The case of an initial coherent state 



When the initial CS evolves under the Hamiltonian of Eq. ( |2.2| ), the (time- 
dependent) expectation values on the RHS of Eq. ( j4.12|) can be evaluated flSSfl 
using Eq. ( p.21|) . The result is the expression for D q (t) given in Eq. ( |C.6| ), 



Appendix C. 

We are interested, in particular, in examining whether fractional revivals are 
accompanied by any significant degree of squeezing and higher-order squeezing. 
For this purpose, we focus on D q (t) at a 2-sub-packet fractional revival. Setting 
t = tt/(2x) in Eq. ( |C.6| ) (Appendix C) and simplifying, we get 

D q {T TeY /2) = ^lyyj 1 + (-l) 9_1 cos 2q9 - 2e- 4l - sin2 ^ 2 cosV}. (4.13) 

As before, 9 is the argument of a (= u 1 / 2 e id ). It is evident that D q vanishes at 
this instant for all even values of q, implying that no even-order squeezing of the 
state accompanies this fractional revival. For odd values of q, however, we find 

D q {T iev /2) = ( sin 2 q6 - e" 4 * cos 2 qO) . (4.14) 

Thus squeezing (or higher-order squeezing) occurs at this instant, provided D q < 
0, i.e., | tan q0\ < e~ 2v . We illustrate this in Figs. [4. l| (a.) and (b), where D\ and 
D 3 are plotted as functions of v for an initial CS with 8 = 0. (Once again, we have 
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Figure 4.1: D q (^T rcv ) versus v for an initial CS (with 9 = 0), for (a) q = 1 and 
(b) q = 3. (Note the different ordinate scales in the two cases.) 

set x = 5 in all the numerical results presented in this Chapter.) We have also 
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Figure 4.2: D q (^T rev ) versus 9 for an initial CS (with v = 0.1), for (a) q 
(b) q = 3. 



1 and 



plotted D\ and D ?J as functions of 9 for a fixed value of v (= 0.1) in Figs. |4~2|(a) 
and (b), showing how squeezing occurs for certain ranges of the argument of a, 
when D q becomes negative. 

An interesting feature is that the CS spreads and loses its coherence almost 
immediately after it starts propagating in the medium. This is borne out by the 
fact that it becomes squeezed even within a very short time of propagation [55| as 
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seen in Fig. [4.3| , which is a plot of D q versus t/T rev for an initial CS with v = 10. 
Further, we see that close to t — the state not only displays squeezing, but also 
higher-order squeezing, for both odd and even integer values of q. 

o.s 

0.6 

0.4 



D 



0.2 

q o 

-0.2 - 

-0.4 - 

-0.6 




1 


1 


/I 

f 

r 

j 

. "r ~ . 
J 

( ' • 




- 


; 


1 j 

J 

/ * 
i i 

i / 


- 


f - 
1 


i 


r- — - 


i 



0.004 0.008 0.012 

t/T 

rev 



Figure 4.3: D q versus t/T rev for an initial CS with v = 10, for q = 1, 2, 3 and 4. 



In the light of these results, we are in a position to investigate the dynamics 
of an initial PACS and compare it with that of an initial CS. 



4.3 The case of an initial PACS 

We show that if the initial state departs even marginally from coherence, as in 
a PACS with a small value of m, the results of the preceding Section change 
significantly. Writing D q as Dq™^ (for ready identification) when the expectation 
values in Eq. (|4.12|) are evaluated for an initial state \a,m), we find the following 
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general result: 



1 m ( 

n=0 ^ 



m + 2q\ v n+q 



n + 2q J n 



x cos 



2(2m + 2n + 2q - l)xqt + v sin Axqt - 2q9 



n + q J n\ 



n=0 

2 



x cos j^(g — 1 + 2m + 2n)xqt + v sin 2%gt — q6 | 

+ E ( q ) , ^ ^ LU-u), (4.15) 
^ \n J (m - q + n)\ mV y v 7 



where n m j n = max (0 , q — m). The derivation of the expression above is given in 
Appendix C. It is the generalization of the expression for D q = obtained for 
a coherent state. As before, we examine at t — |T rev for possible squeezing 
and higher-order squeezing of the state at that instant. The foregoing expression 
reduces at this instant of time to 



2g — 2f(l— cos qff) 



L m (-u 
q 

q\ m 



2&((-l)«i/)} cosV 



71 — T^min 



^ \n/(m-g+n)! 



(4.16) 



We use this to analyze various cases numerically. In contrast to what happens for 
an initial CS, it turns out that there is no odd-power amplitude-squeezing for an 
initial PACS. Even-power amplitude-squeezing does occur, though, for sufficiently 
large values of v. This is illustrated in Figs. |4.4| (a) and (b), which show the ranges 
of v for which falls below zero. Figures |4.5| (a) and (b) depict the variation 
of D% with the phase angle 9 for a fixed value of v. 



Turning to the extent of squeezing as a function of time, in Fig. |4.6| we compare 
the temporal variation of the standard deviation Ax for an initial CS (m = 0) 
and an initial PACS (m = 1), with a = 1 (hence 6 = 0). The horizontal dashed 
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line demarcates the level below which the state is squeezed. It is evident that 
squeezing in the rigorous sense accompanies the fractional revival at t — \T VCV 
when the initial state is a CS. This feature is suppressed when it is a PACS, 
although Ax does dip down considerably around this fractional revival. 

To complete the picture, we have also plotted, in Figs. [17| (a) and (b), Dq 
versus t/T rev for the initial states \a, 1) and \a, 5) respectively, for v = 10 and 
q = 1, 2, 3 and 4. 

We have focused on g th -power amplitude-squeezing at fractional revivals, as 
this turns out to provide rather more discriminatory signatures of higher-order 
squeezing effects than the other alternative, namely, Hong-Mandel squeezing. 
However, a few remarks on the latter are in order here. The relevant variables in 
the case of Hong-Mandel squeezing are 



a + a 



x q and 



a — a 



p q . (4.17) 



For q — 1, of course, Hong-Mandel squeezing is the same as amplitude-squeezing, 
but the two kinds of squeezing differ for q > 2. Figures |4.8| (a,) and (b) show how 
((x — (x)) A ) = ((5x) 4 ), the fourth moment of x about its mean value, varies over a 
revival period for an initial CS and PACS, respectively. The horizontal dotted lines 
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Figure 4.7: .Dg" 1 -' versus t/T rev for an initial PACS (a) \a,l) and (b) \a, 5), with 
f = 10, for g = 1, 2, 3 and 4. 



indicate the bound on ((Sx) 4 ), below which fourth-order Hong-Mandel squeezing 
occurs in this quadrature. An initial CS exhibits such squeezing near revivals, and 
comes close to doing so near the fractional revival at |T rev , but does not actually 
do so. An initial PACS does not display such higher-order squeezing at any time, 
although ((5x) A ) attains its lowest value at revival times. However, fractional 
revivals are marked by rapid oscillations of ((Sx) 4 ), these being most pronounced 
around the 2-sub-packet fractional revival. These features are enhanced further 
in the case of initial states with larger values of m. 



4.4 The Wigner function and the non-classicality 
indicator 

Finally, we examine the Wigner functions corresponding to the wave packets at 
instants of fractional revivals, to quantify non-classical behavior during their time 
evolution. It has been suggested |30| that the "extent" to which the Wigner 
function becomes negative (as a function of its complex argument) may serve as 
a measure of the non-classicality of the state concerned. The normalized Wigner 



function W((3 ; t) (where (3 6 C) is defined in the case at hand as follows fSBfl . One 
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Figure 4.8: ((5x) 4 ) versus time in units of T rev for initial states (a) |a) and (b) 
\a, 5), with v — 1. Note the very different ordinate scales in the two cases. Hong- 
Mandel squeezing in x occurs below the dotted horizontal line. 
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first defines the symmetric characteristic function 



K(z, t) = Tr [p(t) e zat - z * a ] , (4.18) 

where zeC, and p(t) is the density matrix. Then W(/3 ; t) is the two-dimensional 
Fourier transform given by 

W(J3 ; t) = ^ J d 2 z K(z, t) e za "~ z ' a . (4.19) 

Using the representation for the density matrix in the Fock basis, i. e., 

oo 

P(t)= 10 R (4-20) 



l,n=0 



one arrives at the following representation for the Wigner function in terms of the 
density matrix elements pi n (t) fl5T| : 



f-. s OO OO 

x(2/3r" i A n(t)^r'(4|/3| 2 )). (4.21) 



For an initial coherent state |a), we have the standard result 



ct o; u,i2 



p ln (0) = -== e -W . (4.22) 
V(! n! 



The corresponding Wigner distribution is given by the well-known expression p8" 



0) = - e - 2 l a ^l 2 . (4.23) 

7T 

This is positive definite everywhere in the complex /3-plane, justifying the appel- 
lation "classical" for an oscillator coherent state \a). The time- dependent density 
matrix element in this case is given by 

pin (t) = ?-JL e -l°l 2 e -«PP-i)-n(«-i)]t. (4.24) 
vZ! n\ 
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Using this in Eq. ([4.21 ), we compute W(j3,t) numerically. Figures |4.9| (a) and 



(b) show how the Wigner function for an initial CS behaves at the instants of the 
2-sub-packet and 3-sub-packet fractional revivals, respectively. 

For the initial photon-added coherent state \a,m), one finds 



Pln{0) ~ m\L m (-u) (l-m)\(n-m)\ ' (4 ' 25) 
Correspondingly, the Wigner function at t — can be expressed in the closed form 

W((3; 0) = T [ / L m (\2{3-a\ 2 )e- 2 ^ 2 . (4.26) 

It is to be noted that this is no longer positive definite for all complex (3, reflecting 
the fact that this initial state is no longer "totally" classical, as it departs from 
perfect coherence due to the photons that have been "added" to \a) to produce 
the PACS. The time-dependent density matrix element in this case is given by 



Pln (t) = , ° " ^ e -4'(-')-M', (4.27) 

ml L m {—v) [I — my. [n — my. 

where l,n > m. (The matrix element pi n (t) vanishes identically for l,n < m.) 
Once again, using this in Eq. (|4.21| ), we compute the Wigner function numerically. 
Figures |4.10| (a), (b) and (c) are plots of W(/3] t) for the case m = 1 at t = 



0, |T rev and |T rev , respectively. The corresponding plots for the case m = 10 are 
shown in Figs. |4.11| (a) , (b) and (c). With increasing m, the oscillations of the 



Wigner function in the /5-plane between positive and negative values become more 
pronounced, and the region of non-classicality becomes more extensive. 

To get an idea of the degree of non-classicality as a continuously varying func- 
tion of time for each of the different initial states we have considered, it is instruc- 
tive to consider the non-negative quantity S defined as pOf 



5(t) = jd 2 (3(\W((3;t)\-W((3;t)^ 

= Jd 2 ^^^; t)\ - 1. (4.28) 
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Figure 4.9: Wigner function corresponding to an initial state |a) with a — 1, at (a) 
t = |T rcv and (b) t = |T rev . Here, and in the succeeding figures, Pi — Re f3 , @ 2 — 
Imp. 



52 




53 




54 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.S 0.9 1 

t/T 

rev 

Figure 4.12: Non-classicality indicator 5 versus time in units of T rev , for initial 
states \a,m) with a = 1. 



The larger the value of 5, the greater is the extent of non-classicality of the state 
concerned, although of course 5 alone does not give a complete picture of the 



oscillations of the Wigner function. In Fig. [4. 12| we have plotted 5 versus t for 
initial states \a) , \a, 1), and |a, 10) with a — 1. It is clear that in the interval 
between t = and t = T rev , S has the most pronounced local minimum at the 
2-sub-packet fractional revival, followed by the local minima at the 3-sub-packet 
and 4-sub-packet fractional revivals. This feature becomes increasingly prominent 
for larger values of m, showing that the extent of non-classicality also increases 
with m, as expected. 



4.5 Concluding remarks 

We have shown that the crucial difference between an initial CS and an initial 
PACS propagating through a Kerr-like medium arises at the instant |T rev •, when 
the former displays squeezing and the latter does not. A possible experimental 
set-up to examine the squeezing properties of the state at |T rev involves passing 
it through a Kerr medium with a known value of x (the third-order nonlinear sus- 
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ceptibility), and therefore T rev , and using the time evolved state at |T rev as the 
input field. Superimposing this on the field from a local oscillator, and performing 
balanced homodyne detection to measure the variance in x, would determine the 
extent of squeezing. Thus, the successful experimental generation of non-classical 
states of light together with the observation of higher-order moments p9| of ap- 
propriate quadrature variables would enable the probing of the rich structure of 
wave packet dynamics in a nonlinear medium. 
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CHAPTER 5 



Ergodicity properties of entangled two-mode 

states - I 



5.1 Introduction 

In the preceding chapters, we have examined the manner in which signatures of 
revivals and different fractional revivals of a wave packet manifest themselves in 
the dynamics of the expectation values of appropriate operators. We have shown 
that, by tracking the time evolution of various moments of certain operators, se- 



lective identification of different fractional revivals can be achieved |18|. Again, 
since the precise nature of the initial state concerned plays a crucial role in deter- 
mining its subsequent dynamics and the non-classical features it exhibits, these 



signatures help assess the extent of coherence of the initial state fl28|. Further, 
it has been shown that the squeezing and higher-order squeezing properties of 
certain quadrature variables in the neighborhood of a two sub-packet fractional 
revival of the wave packet provide quantifiable measures of the degree of departure 
from coherence of the initial state |29| . 



In this chapter, we extend these investigations to the dynamics of entangled 
states. To be specific, we consider a physically relevant generalization of the 
example considered earlier, and now study the dynamics of a single-mode electro- 
magnetic field interacting with the atoms of a nonlinear medium [34], taking into 



account the interaction between the field mode and the medium. Although we 
consider initial states that are direct product states of the field and atom modes, 
entanglement of these two modes sets in during temporal evolution. However, 
we have shown that for certain ranges of values of the parameters and coupling 
constants in the governing Hamiltonian, the modes disentangle at specific instants 
during the evolution, and the state of the system returns to its initial form, apart 



from an overall phase. This feature is clearly the analog of wave packet revivals 
in the dynamics of a single mode. The entropy of entanglement of the system 
provides a measure of the extent of revival of the original state. At the instants of 
revivals (these being integer multiples of a revival time which we refer to as T rev 
in this two- mode example as well), this entropy vanishes. 

It is evident that, as in the single-mode case, the dynamics of the quantum state 
between t — and t = T rev would lead to closed orbits in a classical phase space 
in which quantum mechanical expectation values are regarded as the dynamical 
variables. The entropy of entanglement is also seen to reduce significantly in 
certain cases at fractions ~, ~, | and | of T rev , signaling the appearance of the 
counterparts of fractional revivals at these instants. We identify suitable operators 
whose expectation values carry clear signatures of this phenomenon. 

We have also investigated the role of the (non-entangled) initial state, and the 
extent of departure of the initial field mode from coherence, on the subsequent 
revival properties of the state. In particular, the link between entanglement of 
states and their squeezing properties |6(| has been examined. The entropy of 



entanglement at any instant increases with increasing departure from coherence 
of the initial field mode, and revivals do not occur. This aspect is studied by 
means of a time-series analysis of the expectation value of a certain observable. 
The details of this analysis and the results obtained are presented in the next 
chapter. 



In the next section, we describe the model |34| we use to study the dynamics of 



two- mode entanglement. Earlier results on the collapse and revival phenomena in 
this case, and their dependence on the initial state and parameter values, are sum- 
marized. In Section 5.3, we show how the extent to which a state revives at any 
instant is related to the instantaneous sub-system von Neumann entropy (SVNE) 
and the sub-system linear entropy (SLE). Certain operators whose expectation val- 
ues carry signatures of the different fractional revivals are also identified. We have 
examined the dynamics of several initial states, taking the oscillator representing 
an atom of the medium to be in the ground state while the field is, respectively, 
in a Fock state, a CS, and a PACS. This enables us to analyze systematically the 
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effects of different initial field modes on the dynamics. The relationship between 
the squeezing exhibited by the state of the system and the initial field mode is 
also brought out. We conclude the chapter with some brief comments. 



5.2 Single-mode field in a nonlinear medium: 
The model 



Sanz et al. [B3 have examined the non-classical effects that arise in the dynamics 



of two entangled modes governed by a nonlinear Hamiltonian, in the framework 
of an exactly solvable case: two modes of an electromagnetic field interacting in 
a Kerr- like medium. Taking the initial state to be a direct product, either of two 
Fock states or of two coherent states, the periodic exact revival of these states has 
been established, and the manner in which these properties show up in the collapse 
and revival phenomena displayed by the expectation values of observables has been 
investigated. The collapses are marked by expectation values remaining constant 
over a certain interval, while the revivals are signaled by rapid pulsed variations of 
these quantities. This is clearly analogous to the signatures of revivals displayed 
in the single-mode case that we have discussed in the preceding chapters. 

As stated in Chapter 1, our objective is to investigate the full range of regular 
and irregular dynamical behavior exhibited by quantum mechanical expectation 
values. For this purpose, we require a system where revivals, near revivals, as well 
as no revivals of the initial state can all occur, for different parameter values. As 
already mentioned, a good candidate Hamiltonian for our purposes is the one that 
describes the interaction of a single-mode field of frequency to with the atoms of 



the nonlinear medium through which it propagates. The medium is modeled [34 
by an anharmonic oscillator with frequency ujq and anharmonicity parameter 7. 
The Hamiltonian of the total system is given by 

+ 7 6 t V + ^(a t 6 + 6 t a). (5.1) 

a and (v are the annihilation and creation operators pertaining to the field, while 
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b and are the corresponding atomic oscillator operators. The coupling constant 
g is a measure of the strength of the coupling between the field modes and the 
atom modes. The Fock basis is given by {\n) a <S> \n') b }, where n and n' are the 
eigenvalues of a) a and tfb, respectively. The photon number operator N = a) a no 
longer commutes with H. However, the total number operator N to t = o)a + tfb 
does commute with H. Hence we may write the basis states as \N — n) a ® \n) b , 
using iV to label the eigenvalue of N tot . (We remark that, in the discussion of 
the one-mode case in the preceding chapters, we have used N in places to denote 
the eigenvalue of the photon number operator N. No confusion should arise as 
a result, as we specify the notation explicitly wherever required.) For notational 
simplicity, let us write 

\N - n) a <g> \n) b = \N - n; n) . (5.2) 

It is evident that (N — n; n\ H \N' — n' ; n') = 0, if iV ^ N'. Hence, for each 
given value of N, the Hamiltonian H can be diagonalized in the space of the 
states \N — n; n), where n — 0, 1, . . . , N. Let the eigenvalues and eigenstates 
of H be \n s and \ipNs), respectively, where s — 0, 1, . . . , N for a given N, and 
N — 0, 1, . . . ad inf . It is convenient to expand \ipNs) in the basis {\N — n; n)} 
as 

N 

\^Ns)=Y. d n S \N-n;n), (5.3) 

n=0 

so that d^ s = (N — n; n\ tpN S ) . 

The time-dependent density operator p(t) corresponding to different choices of 
the initial state ^(O)) can now be obtained. For instance, in the simple case in 
which |^(0)) = |JV) a ® |0) 6 = \N; 0), we have 

N N 

p(t) = exp [-i(\ Ns - \ Ns >)t] (ip Ns \N ; 0) 

s=0 s'=0 

x (N ; 0\ip Na >) \iiNs) {ipNs'\ ■ (5.4) 

Such expressions for the time- dependent density matrix and the corresponding 
reduced density matrices for the two sub-systems (i. e., the field and the medium, 
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respectively) are required in order to evaluate the expectation values and entropies 
to be considered in the next section. Further details are given in Appendix D. 



We now summarize some of the key results of Ref. [34| pertaining to the 
collapse and revival phenomena exhibited by the field, as these are of relevance 
to our study When g = 0, the Hamiltonian in Eq. ( |5.1| ) is in diagonal form as 
it stands. If 7 = 0, H is exactly diagonalizable, and there is periodic exchange of 
energy between the field and the atomic oscillator. For non-zero values of the ratio 
7/(7 of the respective strengths of the nonlinearity and the field-atom interaction, 
collapses and revivals could occur over certain intervals of time, in between these 
periodic exchanges of energy. This phenomenon translates into the behavior of 
expectation values of certain observables as well. For instance, during a collapse 
of the energy exchange over an interval of time, the mean photon number (a^a) 
essentially remains constant. A revival of the energy exchange is signaled by rapid 
oscillations of the mean photon number about this steady value, over the relevant 
time interval. 

Further, when the atomic oscillator is initially in its ground state, while the 
field starts either in a Fock state or in a coherent state, one finds the the following 
results: (a) For weak nonlinearity (7/5 <C 1), collapses and revivals of the mean 
photon number occur almost periodically in time, for both kinds of initial field 
states. The revival time can be shown to be approximately equal to 27r/7 in the 
former case, and 47r/7 in the latter, (b) For j/g ~ 1, such collapses and revivals 
occur more irregularly if the field is initially in a coherent state, compared to the 
case when it is initially in a Fock state. (c) For large nonlinearity (j/g 1), 
these collapses and revivals do not occur. Bearing these results in mind, in the 
next section we examine the manner in which the foregoing collapse and revival 
phenomena are mirrored in the entropy of entanglement of the system. We identify 
observables which carry signatures of these collapses and revivals, and discuss the 
influence of the departure from coherence of the initial state of the field on the 
extent to which it revives subsequently. 
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5.3 Entanglement properties 



Let us now examine the detailed dynamics of three different classes of initial 
states which are direct products of the field and atom states, evolving under the 
Hamiltonian in Eq. ( |5.1|) . As stated earlier, the initial state of the atom is taken 
to be the oscillator ground state |0) 6 , while that of the field is, respectively, (a) 
a Fock state \n) a , n — 0, 1, . . . ; (b) a CS \a) a ; and (c) an m-photon-added 
CS \a,m) a , m — 1, 2, ... . (Recall that the suffixes a and b correspond to the 
electromagnetic field and the atoms of the medium, respectively.) Extending the 
notation in Eq. ( |5.2|) , we shall write 

\a,m) a ® \n) b = |(a, m) ; n) (5.5) 

for the direct product state concerned. 

As we are dealing with a pure bipartite system, it is natural to consider the 
time-dependences of S k , the sub-system von Neumann entropy (SVNE), and 5 k , 
the sub-system linear entropy (SLE), where the suffix k stands for either a or b, 
depending on the sub-system considered. These quantities are defined as 

S h (t) = -Ti k (p k {t) In p k (t)) (5.6) 

and 

5 k (t) = l-Tr k (pl(t)), (5.7) 

where p k (t) is the time-dependent reduced density operator for the sub-system 
concerned. In terms of the set of eigenvalues of p k (t), we have 

S k (t) = ^ Ai i} (t), 5 h (t) = [A«(t)] 2 , (5.8) 

i i 

where the summation runs over all the eigenvalues X h (t). Once again, as we have 
a pure bipartite system, the SVNEs for the two sub-systems are in fact equal to 
each other at all times, as are the SLEs. (This used as one of the checks on the 
numerical routines employed.) In Appendix D, we have outlined the steps leading 
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to the derivation of the matrix elements of the reduced density matrices p a {t) and 
Pb(t), required for the numerical computation of Sk(t) and Sk(t). 

We are interested, in this chapter, in examining revival phenomena. We there- 
fore restrict ourselves to the case of weak nonlinearity, i. e., j/g << 1, as this 
is the situation in which these phenomena occur. For illustrative purposes, we 
set the values of the parameters at u — ujq — 1, and 7 = 1, g = 100, so that 
j/g — 10 -2 , in all the numerical results reported in the rest this section. Note 
that the revival time then corresponds to gt ~ 2007T for an initial Fock state of 
the field, and gt m 400n for an initial coherent state of the field. 

Figures |5.1| (a) and (b) depict plots of the SVNE and the SLE versus gt for 
respective initial states |10 ; 0) and \a) a <8> |0) 6 = \a ; 0) with the parameter value 
v = 1. The band-like appearance of the plots in Fig. 5T (a) arises from the 
extremely rapid oscillations of the ordinates. The corresponding plots for an initial 
state \a, 5) a <8> |0) 6 = \(a>, 5) ; 0) in the cases v — 1 and v — 5 are shown in Figs. 
|5.2| (a) and (b) respectively. In all the cases above, the SVNE (the upper plot in 
each figure) is larger than the SLE (the lower plot in each figure) at any instant 
of time. It is evident that both SVNE and SLE display roughly similar oscillatory 
behavior in time. However, certain striking differences arise in the time evolution 
of the SVNE and SLE, depending on the actual initial state considered. If the field 
is initially in a Fock state or a CS, the entropies return to zero at regular intervals 
of time (see Fig. |5.1| ), signaling a revival of the initial state. (Note that the revival 
times are indeed approximately equal to 2n and 4ir, respectively, recalling that 
we have set 7 = 1 and g = 100.) In contrast, if the initial state of the field is 
a PACS, the entropies do not exhibit such exact revivals, although they do show 
substantial oscillatory behavior near T rev . Further, with an increase in the value 
of is, the oscillations in the SVNE and SLE die down. This effect is enhanced for 
larger values of m, as seen in the rapid increase and saturation of both the SVNE 
and SLE for an initial state \(a, 5) ; 0) for v — 5, in contrast to the corresponding 
plots for is — 1: see Figs. [5]2] (a) and (b). 



It is also clear that the SVNE and SLE display marked oscillatory behavior 
near |T rev , |T rev and \T rev . This behavior may be regarded as the counterpart, 
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Figure 5.1: SVNE and SLE versus gt for (a) an initial Fock state |10; 0), and 



(b) an initial coherent state \a ; 0) with v = 1. The ratio j/g = 10 2 in Figs. O 
through |5.11. 
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in the coupled system, of the fractional revivals seen in the single-mode case (as 
described in Chapters 2 and 3). Again, these oscillations die down in amplitude 
with increasing m when the initial state of the field is a PACS, and are most 
pronounced if the field is initially in a Fock state. 

We now examine whether signatures of the features seen above appear in the 
time evolution of expectation values of observables. For this purpose, we define 
the quadratures 

£= (x a + x b )/2, r 1 ={ Pa + Pb )/2, (5.9) 

where 

x a = (a + atyVH, x b = (b + b ] )/V2 (5.10) 

and 

Pa = (a - J)/{iV2), p b =(b- tf)/{iV2). (5.11) 

If the field is initially in a Fock state, both (£) and (rj) vanish identically at all 
times. Setting u — ujq — 1, 7 = 1 and g = 100 as before (for ready comparison 
with the behavior of the SVNE and SLE discussed above,), we have plotted (£) 
versus gt for an initial state \a ; 0) with v — 1 and 5, respectively, in Figs. (a) 
and (b). Figure |5]3| (a) shows that (£) displays rapid pulsed oscillations near t = 
Att, similar to its behavior near t = 0. This mark of a revival is consistent with the 
behavior of the SVNE and SLE in this case, cf. Fig. |5.1| (b). The collapses are not 
sharp, in the sense that (£) is not constant over the time interval between successive 
revivals — oscillatory bursts occur in between, with a slight enhancement of these 
oscillations around the fractional revival at \T rew . In contrast, the collapses in 
between revivals are much more complete for larger values of u, as seen in Fig. 
|5.3| (b), consistent with the corresponding behavior of the SVNE and SLE in this 
case. 

An interesting feature is that these collapses become much sharper for even 
a marginal departure from coherence of the initial state of the field. (£) remains 
virtually constant over the duration of the collapse, and then bursts into rapid 
oscillations close to revivals, as seen in Fig. [5.4| which corresponds to m = 1 
and v — 1. As in the case of single- mode dynamics, the amplitude of the oscil- 
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Figure 5.4: (£) versus gt for an initial state |(a, 1) ; 0) with v — 1. 

lations in the neighborhood of T rev decreases significantly with an increase in m. 
Thus, for small values of u, it is easy to distinguish between an initial CS and 
an initial PACS. The expectation value of rj also displays these signatures. While 
the sub-system variables x a , Xb, p a and pb do exhibit revivals in their expectation 
values, their higher moments do not capture the occurrence of fractional revivals. 
However, the higher moments of the combinations £ and r\ do carry distinguishing 
signatures to selectively pin-point the analogs of the different fractional revivals 
that occur in the single-mode case. Hence £ and 77 are the appropriate dynamical 
variables for the interacting system under consideration. 

The standard deviation Z\£ of £ reflects the occurrence of the dips in the plots 
of the SVNE and SLE at \T rev . Consider, first, the case when the initial state is of 
the form \N; 0). Then a plot of Z\£ versus gt shows a burst of rapid oscillations at 
t w 7T (pa |T rev ). Fig. U illustrates this for N = 10. (Note that the case iV = 1 is 
exceptional in this regard, since A£ reduces to a constant in this instance.) This 
feature holds for an initial CS or PACS as well, as is evident (see Fig. [5.6| ) from 
the sudden burst of oscillations in Z\£ around t ~ 2n ps |T rev for initial states 
\a ; 0), I (a, 1) ; 0) and \(a, 5) ; 0) (recall that T rev pa 4n in this case). 

We note that squeezing occurs in the neighborhood of |T rev when the initial 
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Figure 5.5: A£ versus gt for an initial state 1 10 ; 0). 




Figure 5.6: Z\£ versus gt for an initial state \(a,m) ; 0) with v = 1 and m — 0,1 
and 5, respectively. 
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state of the field is a coherent state: drops below the value \ (the horizontal 



dotted line in Fig. |5.6| ) in the case when the initial state is |ot; 0), in contrast 
to what happens for an initial PACS |(a, m); 0). While this is similar to the 
squeezing property seen in the single-mode case, such parallels do not hold in 
the case of higher-order squeezing. The relevant quadrature variables in this case 
(analogous to those in Eq. (|4.7|) ) are 

_ (a 9 + at* + \fl + ftt9) ( a ? _ a t? + ftfl _ &t?) 

Zl " 27! ' 2 " Wf ' ( } 

For q = 1, Z\ and Z 2 reduce to £ and 77, respectively. The indicator of higher-order 
squeezing is defined, analogous to D q (t) in Eq. (|4.11| ) of Chapter 4, as 

(AZ,f-\\(\Z 1 ,Z 2 })\ 

DM = Ij z 2 ]> I • (5 ' 13) 



where denotes the standard deviation in the variable Zi defined in Eq. ( p. 12 ). 



We note that \Z\ , Z 2 ] is no longer a polynomial function of the (total) number 
operator N tot alone. Extending the notation used in Chapter 4, we shall write 
D { q m) for D q in the case of an initial state \(a,m) ; 0), and D q N ^ to denote D q 
for an initial state \N ; 0). In contrast to the single-mode case, even for weak 
nonlinearity (j/g = 10~ 2 ) and v — 1, amplitude-squared squeezing (q = 2) is 
absent at t — |T rev whether the field is in a Fock state, or a CS, or a PACS. For 
iV = 10, and an initial state |10; 0), we have plotted D2 versus time in Fig. 

The state is never squeezed, as D2 is always positive. Similar plots for an 
initial state \a ; 0) and \(a, 1) ; 0) are given in Fig. |3 While D ( 2 0) is almost zero 
at certain instants for a coherent initial field state \a) a , it is clear that it does not 
drop below the horizontal line indicating the value below which amplitude-squared 
squeezing occurs. 

Turning to the higher moments of £ and r], we note that all the odd moments of 
£ vanish identically for all t if the initial state is a direct product of Fock states (see 
Appendix D). Analogous to the single-mode case, for small values of v, the higher 
moments of £ show distinct signatures at fractional revivals only if m is sufficiently 
large. This is evident from Fig. p^9| , where we have plotted ((£ — (£)) 3 ) = ((<5£) 3 ) 
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Figure 5.7: versus gt for an initial state \N ; 0) with N = 
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Figure 5.9: ({SC,) 3 ) versus gt for an initial state (a) \a ; 0) (b) \(a, 3) ; 0) with v — 1 
and j/g = 10~ 2 . 



versus time for initial states \a; 0) and \(a , 3) ; 0), respectively, with v = 1 and 
j/g = 10 -2 . In the latter case, the signatures are manifested in the pronounced 
dips in ((<5£) 3 ) at |T rev and |T rev . However, for larger values of u, such signatures 
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Figure 5.10: ((<5£) 3 ) versus gt for an initial state \a ; 0) with v = 5. 



appear even in the case of an initial CS (m = 0), as is evident from Fig. |5.10| . In 
contrast to this, the variances and higher moments of the subsystem variables x a , 
Xb, p a and pb do not display these signatures. For instance, the plot of Ax a versus 
time shown in Fig. |5.11| does not display a distinctive signature at |T rev alone: 
apart from oscillations in the vicinity of this instant, we see marked dips around 
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|T rcv and |T rev as well. 
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Figure 5.11: Ax a versus gt for an initial state \a ; 0) with v — 5. 

5.4 Concluding remarks 

In the light of our discussion, it is clear that the quadrature variables £ and rj 
corresponding to the full system are the appropriate observables whose dynamics 
signals the occurrence of the analogs of revival phenomena, while the expectation 
values of sub-system variables would seem to be less suitable choices for high- 
lighting collapse and revival phenomena. On heuristic grounds, we may, however, 
expect them to be better candidates for investigating non-periodic temporal be- 
havior — in particular, time evolution that is actually indicative of exponential 
instability. We address this issue in the next chapter, in which we analyze the dy- 
namics of the sub-system variable (a^a) for a wide range of values of the parameter 
v and the ratio j/g. 
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CHAPTER 6 



Ergodicity properties of entangled two-mode 

states - II 



6.1 Properties of the coupled-system Hamilto- 
nian 

In this chapter, we continue our investigation of the dynamics of a single-mode 
electromagnetic field propagating in a nonlinear medium, with the model Hamil- 
tonian of Eq. ( |5.1|) . We write down H again for ready reference, restoring the 
over-all factor of h for the moment, as in Eq. (|1.9|) : 



+ 7 &t 2 & 2 + 0(o t 6 + & + a)]. (6.1) 

We recall that the operators a and a) pertain to the sub-system represented by 
the field, while h and are the ladder operators for the other sub-system, the an- 
harmonic oscillator modeling an atom of the medium. We have found in Chapter 
5 that sub-system variables do not adequately capture the effects of revival phe- 
nomena. On the other hand, this would make them potentially good candidates 
for examining the ergodicity properties of the system. The mean photon number 
(N) = (aJa), which is also essentially the average energy of the field mode, turns 
out to be an ideal candidate for this purpose. Our reasons for this choice are as 
follows. 

In the single-mode model with Hamiltonian K% N(N — 1) considered in Chap- 
ters 2, 3 and 4, it is evident that the expectation value of N, as well as all its 
higher moments, are constants of the motion for any initial state. In contrast to 



this, N does not commute with the Hamiltonian ( |6.1|) of the coupled system. As 



a consequence of the entanglement between the atom and field modes, the mean 



photon number is no longer a constant of the motion for any non-zero value of 
the coupling constant g. However, as already stated, the total number operator 
Ntot = + b'b commutes with H for all values of the parameters in H. Its ex- 
pectation value is therefore always a constant of the motion. In the absence of the 
anharmonicity parameter 7, the coupled Hamiltonian is essentially linear in each 
of the sub-system variables, and can be diagonalized in terms of linear combina- 
tions of the original ladder operators. In physical terms, this leads to fairly simple 
dynamics, essentially entailing a simple periodic exchange of energy between the 
two sub-systems or modes. (This mechanism becomes even more transparent if we 
set u) — Uq .) When 7 is non-zero, the dynamics is more complicated, but the ex- 
istence of the operator N tot ensures that the system as a whole is "well-behaved" , 
with a discrete spectrum (labeled, as we have seen, by the quantum numbers iV 
and s, where N = 0, 1, . . . and s = 0, 1, ... , N). The time- dependence of (N) is 
thus a direct consequence of the coupling of the field mode to another degree of 
freedom; and the deviation of (N) from periodic temporal variation is a probe of 
the effects of the nonlinearity in the latter. It turns out that a remarkable diversity 
of temporal behavior is exhibited by (N), in different ranges of the parameters in 



In this context, it is very helpful to understand the dynamical behavior of the 



also needed in order to facilitate subsequent comparison of the quantum dynamics 
with the classical one. We may obtain the classical Hamiltonian as follows. Let 
the linear harmonic oscillator associated with a and a' have a mass m, position 
x and momentum p x . Similarly, let the oscillator associated with b and have a 
mass M, position y and momentum p y . Putting in all constant factors (including 
h), we have 



The hermitian conjugates a' and are of course obtained by replacing % with 
— i in the above. Substituting the foregoing in Eq. (|6.1|) , we get H in terms 
of x, p x , y and p y . Passing to the classical limit in which h — > 0, we find that 



H. 



completely classical counterpart of the Hamiltonian of Eq. ( |6.1| ) . This analysis is 




(6.2) 
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the only consistent way to obtain a non-trivial, finite expression for the classical 
Hamiltonian in this limit is to let 7 also tend to zero simultaneously, such that 
the ratio j/H tends to a finite value A, say. When this is done, the classical 
Hamiltonian obtained is 



2m 2 2M 2 ° y u% \ 2M 2 



,2„,2 



+ -^^VrnM^o^+^LV (6.3) 

with the canonical Poisson bracket relations {x , p x } = {y , p y } = 1. As it stands, 
it is not immediately obvious that the corresponding motion is bounded for all 
initial conditions, owing to the presence of the cross terms xy and p x p y that can 
assume either sign. However, the counterpart of N tot = a^a + tfb, namely, 

Px 1 I. , 1 { Py , 1 a f ,2, ,2 



jV d = _ i2L + -mwV + — — + -^o¥ , (6.4) 
a; 2 J u \2M 2 ° y J ' v ; 

is a second, analytic, constant of the motion ({N c i , H d } = 0). This ensures, of 
course, that the two-freedom system is Liouville- Arnold integrable. Further, since 
N c \ = constant is evidently a bounded hypersurface (an ellipsoid) in the four- 
dimensional phase space, it is clear that motion under H c \ is indeed bounded for 
any set of initial conditions. All four Lyapunov exponents vanish, and the classical 
motion is always regular (and restricted to a 2-torus) for each set of initial condi- 
tions. The dynamics can be analyzed in further detail, including the identification 
of the location and nature of the corresponding critical or equilibrium points. Sev- 
eral features emerge that are interesting in their own right, but we do not go into 
these here, as they are not directly germane to our present purposes. We note, 
though, that the equation of motion of the classical counterpart of N (= a^a), 
on evaluating the Poisson bracket concerned, is 



dJ 



1 f pI , 1 jj 



H — mu x 
lu \ 2m 2 




x Pv ~ V VP* ■ ( 6 - 5 ) 

V mu ' 



The equation of motion of the combination on the right-hand side of Eq. ( |6.5|) 
involves other functions of the dynamical variables, and so on. The resulting set 
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of equations does not close. This is brought out more clearly in the quantum 
mechanical problem, to which we now revert. 

Returning to the quantum mechanical Hamiltonian ( |6.1|) , some additional in- 
sight into its precise nature is provided by re-writing it in terms of angular mo- 
mentum operators ||34|| . The latter are defined in the standard manner in terms 
of the two mutually commuting sets of boson operators, according to 

J+ = a f 6, J_ = a6 f , J z = \(a)a - tfb). (6.6) 



The Hamiltonian in Eq. (|6.1|) becomes 



H = h 



~(u + u - 7 )N tot + (u - uo + i)J z + ^(Ntot - 2J Z ) 2 + 2gJ x 



• (6-7) 



We also have the well known relation J 2 = |N tot (^N tot + l). Therefore, as J 2 
is a function of N tot , it is evident that N to t commutes not only with H, but 
also with every component Jj of the angular momentum. However, [J z , J x ] ^ 0, 
of course, and this suffices to preclude exact diagonalization of H. It is clear 
that we have a problem similar to that of an Ising spin in a transverse magnetic 
field. The quantity whose time evolution we shall track is the expectation value 
of N = a^a = J z + | N to t • Using the standard algebra of angular momentum 
generators, the equation of motion of this operator works out to 

= W = {ihyl [N ' H] = ~ ig {aib " ah]) = 2gJy - (6,8) 



But we also find 



and 



^ = (u - uj - 7 N to t) Jx -2gJ z + -f [J z , J x ] + (6.9) 



(7 Ntot - CO + UJ - 7) Jy - 7 [Jy , J z } + , (6.10) 



dt 

where [• • • , •••]+ stands for the anti-commutator of the operators concerned. 
When differentiated with respect to t, these terms lead to trilinear combinations of 
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the angular momentum operators, and so on. It is easily checked that the system 
of equations does not close. 

Having obtained some idea of the dynamics we can expect of the operator N, 
we now turn to a numerical study. We shall regard the mean photon number (N) 
(equivalently, the mean energy of the field mode) as an effective classical dynamical 
variable, and analyze its time series over long intervals of time. We carry out this 
analysis for a wide range of values of the ratio j/g, which is a measure of the 
relative strengths of the nonlinearity and the interaction present in the coupled 
system. As in Chapter 5, we shall take the atomic oscillator to be initially in the 
ground state, and the field to be in a CS or a PACS. 

In order to set the framework, in the next section we briefly outline the rele- 
vant features of time-series analysis [pql , as this topic is quite distinct from those 
involved in the preceding parts of this thesis. 

6.2 Time-series analysis in brief 

We consider a typical situation, in which the time-series data of just one observ- 
able (signal), or that of a certain sub-set of observables corresponding to a given 
physical system, are available for analysis. In general, there will be other dynam- 
ical variables whose time series are not readily accessible, but which also control 
the full dynamics of the system. The dynamical equations are typically a set of 
nonlinear coupled evolution equations for the complete set of dynamical variables. 
In time-series analysis one attempts to ascertain, from the available data, the ef- 
fective number of variables that are actually relevant in determining the temporal 
behavior of the system, consequently establishing the minimum number of effective 
dimensions (i em b (the embedding dimension) of the "phase space" of these relevant 
variables. Evidently, these variables have to be generated from the values of the 
observed variables at a given time and those corresponding to their time-delayed 
copies. Hence this space is only formally equivalent to the actual phase space 
of the system, and the dynamical variables generated from the data are merely 
effective variables. However, if the analysis is carried out appropriately, they will 
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(in most cases) capture the salient features of the dynamics in the actual phase 
space of the system — in particular, of dynamical chaos, if any. 

The embedding theorem indicates how d em b can be obtained. To begin 



with, one assumes different integer values d for the dimensionality of the effective 
phase space, and generates the d dynamical variables from the time-series data 
available. Each of these variables is regarded as a component of a d- dimensional 
vector. The embedding dimension will turn out to be one of the values assumed 
for d, obtained by the procedure described below. 

Generating the set of effective dynamical variables at any instant involves the 
identification of a time delay r such that, from the available time-series of the 
signal, data points separated by time intervals > r are sufficiently independent of 
each other to be treated as the independent dynamical variables. We assume that 
the time series is obtained by sampling the variable(s) concerned at time steps St. 
Clearly, if two data points s(n) and s(n + T) at discrete times n and (n + T) (in 
units of St) are sufficiently independent of each other, the information deducible 
regarding the measured value of s(n) through the measurement of s(n + T) must 
tend to zero, for sufficiently large T. More precisely, if p(s(n)) and p(s(n + T)) 
are the individual probability densities for obtaining the values s(n) and s(n + T) 
at times n and (n + T), respectively, and p(s(n) , s(n + T)) is the corresponding 
joint probability density, the average mutual information 

KT)= yj ^(").^+r))iog 2 [ f; ( "'-; ( " +J 2| ] (6.H) 

.W^+D Lp(»(n))p(»(n + T))J 

should tend to zero for sufficiently large T. In practice, a frequently-used prescrip- 



tion [37| is to take r to be that value of T at which the first minimum in I(T) 
occurs. Thus, from a time series s(0), s(l), s(2), ... of the signal measured at dis- 
crete times t = 0, 1, 2, one constructs the phase space vector s at t = 0, with 
components s(0), s(r), s(2r), . . . , s((d— l)r). After the first time step, the vector 
s evolves to Si , with components s(l), s(l + r), s(l + 2r), . . . , s(l + (d — 1)t), 
and so on. This procedure is used to reconstruct sets of time-delayed vectors in 
the (i-dimensional phase space. Each vector defines a point in the space. Taken 
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in the correct sequence, these points describe the time evolution of the underlying 
dynamical system. 

This procedure is carried out for different values of d. The next task is to 
identify the value of d that corresponds to <i em b . For this purpose, we evaluate 
the correlation integral 

C(r)= lim — V0(r-|si-s,-|)= / dVc(r'). (6.12) 

i,j=0 J0 

C(r) is an estimate of the average correlation between the various points in a 
phase space of a given dimensionality d that contains a large number of points n 
generated from the time series data. The integrand c(r) itself is the standard cor- 
relation function. It has been argued [36| that, for sufficiently small r, C(r) scales 



like a power of r ; moreover, once d reaches the correct minimum embedding 
dimension d em b , the exponent ( will not change with an increase in d beyond d em b- 
This determines the exact value of d emb . In practice, d em ^ could be much smaller 
than the dimension of the actual phase space of the system considered. However, 
in a generic situation, the dynamics of the effective variables alone, in the phase 
space of dimensionality d em b , should suffice to extract quantitative information 
about the dynamical chaos in the signal. 

A specific objective of our study is to determine whether a sub-system variable 
( an appropriately chosen expectation value) can exhibit exponential sensitivity 
and instability, i. e., "chaotic" behavior. For this purpose we need to calculate 
the maximal Lyapunov exponent in the reconstructed phase space of dimension- 
ality d em b . There are, in general, d em b Lyapunov exponents corresponding to the 
dynamics in the reconstructed phase space. For our purposes, it suffices to esti- 
mate the largest Lyapunov exponent A max of this set. If this quantity turns out 
to be a positive number, we have chaotic behavior. We use a robust algorithm 
|62|1 developed for the estimation of the maximal Lyapunov exponent from data 



sets represented by time series (see also |63|). We start with the set of "distances" 
{dj(0)}, where dj(0) is the separation between the j th pair of nearest neighbors in 
the phase space. This evolves under the dynamics to the set of distances {dj(k)} 
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after k time steps St, i. e., at time t — k8t. The quantity (In dj(k)), where the an- 
gular brackets denote the average over all values of j, is plotted against the time, 
k 5t. The result is, in general, a curve that has a short initial transient region, 
a long, clear-cut linear region, and a subsequent saturation region. The slope of 
the linear region, identified by a least squares fit, is precisely the largest Lyapunov 
exponent sought, A max . Once again, it can be established that, if the procedure for 
phase-space reconstruction is implemented carefully and the minimum embedding 
dimension <i cm b obtained correctly, any further increase in the dimensionality of 
the reconstructed phase space should not alter the inferences made regarding the 
exponential instability, if any, of the system. 

With this summary of the relevant aspects of time series analysis, we proceed 
to carry out a detailed time series analysis of the values of the mean photon 
number obtained at regular intervals of time, as the single-mode field propagates 
in the nonlinear medium. The corresponding power spectra of the signal have 
been obtained using a standard package (Labview). 

6.3 Ergodicity properties of the mean photon 
number 

As already stated, we take the atomic oscillator in its ground state and the field to 
be either in a CS or a PACS, initially. Again, as in Chapter 5, we set uo = uj = 1 for 
illustrative purposes. We first consider the case of weak nonlinearity, i. e., j/g <C 1, 
with v — 1. We recall from Chapter 5 that, if j/g <C 1 and v is sufficiently 
small, then the state exhibits revivals (or approximate revivals) if the field mode is 
initially a CS or a PACS. Distinctive signatures of such revivals show up both in the 
entropy of entanglement and in the dynamics of appropriately chosen quadrature 
variables. At revivals (or near revivals), the expectation values of all observables 
return to their initial values (or to the neighborhood of those values). In a classical 
phase space spanned by these expectation values, the trajectories may therefore 
be expected to be periodic or ergodic (as in the case of quasi-periodicity), but 
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not chaotic, in any case. However, while this indicates that the system as a 
whole (comprising, in our case, the field and atom modes in interaction with each 
other) is a non-chaotic dynamical system, there is no a priori reason to assume 
that this remains true for the field or atom sub-systems taken individually - 
particularly because the sub-system variables do not necessarily capture signatures 
of the revival phenomena, as pointed out in the preceding chapter. Keeping these 
aspects in mind, we have analyzed the time series data of the mean photon number 
N, and examined the plots of the corresponding power spectrum S(f) (which 
is essentially the Fourier transform of the autocorrelation function of the mean 
photon number, / denoting the frequency), in order to make deductions regarding 
the nature of the dynamics. The frequency ratio f/g is clearly the counterpart of 
the dimensionless time variable gt in terms of which the dynamics of the field mode 
has been investigated in Chapter 5. As a representative case of weak nonlinearity, 
we have chosen (as before) the parameter values 7 = 1 and g = 100. The time step 
we have used in the numerical work related to this case is 5t = 10 -2 . Figure |Q| (a), 



which is a plot of the logarithm of the power spectrum versus f/g when the initial 
state is \a ; 0), corresponds to quasi-periodic behavior. With increasing departure 
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Figure 6.1: Power spectrum of the mean photon number for the initial states (a) 
| a ; 0) and (b) \(a, 5) ; 0) with 7/g = 10~ 2 and v — 1. 



from coherence of the initial field mode, the number of frequencies seen in the 
power spectrum increases, as is evident from Fig. |6j] (b), which corresponds 
to an initial state |(a,5) ; 0). Further, for a given initial state, the number of 
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characteristic frequencies in the power spectrum increases with an increase in v. 



This is evident when we compare Figs. |6.1| (b) and Fig. |6.2| , which depicts the 
I 1 1 1 1 1 1 1 1 1 — 
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Figure 6.2: Power spectrum of (N) for an initial state \(a, 5) ; 0) with 7/g = 10 2 
and v = 5. 



power spectrum for the same situation, but with v = 5. The embedding dimension 
turns out to be c4mb = 3 in all these cases, and a detailed time-series analysis 
confirms that the maximum Lyapunov exponent is zero. 

In contrast to the case of weak nonlinearity, the nature of the dynamics of the 
sub-system changes drastically when 7/g > 1. It ranges from quasi-periodicity 
to fully-developed chaos, depending on the precise nature of the initial state. As 
representative values for this nonlinearity-dominated regime, we have set 7 = 
5 , g = 1. We use a time step 5t = 10" 1 in this instance, so that gSt = 10 _1 , as 
this suffices to pick up the oscillations in the variables concerned. We first examine 
the case corresponding to v — 1. For an initial field mode which is a CS, both the 
time series and the power spectrum confirm that the sub-system dynamics is not 
chaotic. In contrast to this, an initial field state that is a PACS leads to chaotic 
behavior, for sufficiently large values of m. The time series analysis yields a value 
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d e mb = 5 in this case. Figure |T3] shows the variation of (\ndj(k)) with t for an 
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Figure 6.3: (In dj(k)) versus t = k5t for an initial state |(a, 5) ; 0) with j/g = 5 
and v = 1. The solid line corresponds to d crn ^ = 5. The dotted lines correspond 
to values of d em b from 6 to 10. 



initial state |(a, 5) ; 0). It is evident that the slope is positive (the actual value 
being w 0.5), indicating that the mean energy of the field mode varies chaotically 
in time. 

In each case, we have increased the dimensions of the phase space beyond 
the embedding dimension, and verified that the value of the maximal Lyapunov 
exponent is a constant independent of this increase. Further, we have verified in 
each case that, if {b^bj is chosen as the signal for which the time series data is 
computed, the behavior is similar to that inferred from the dynamics of (a^a), in 
such a manner that (a^a) + (b^b) is a constant, as required. This lends support to 
the conclusion that round-off or truncation errors may be ruled out as the source 
of the computed chaotic behavior. 

For a given value of 7/(7, an increase in either v or m in the initial PACS 
results in an increase in the value of A max . This is evident from Figs. |0] and 
|6.5| , where plots of (lndj(k)) versus t are shown for initial states 1) ; 0) and 
I (a, 5) ; 0) for v = 5. We find that A max ~ 0.58 and 0.85, respectively, in these two 
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Figure 6.4: (In dj(k)) versus t for an initial state | (ck, 1) ; 0) with ^/g — 5 and 
z/ = 5. The solid line corresponds to embedding dimension equal to 5. The dotted 
lines correspond to values of d em \, from 6 to 10. 
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Figure 6.5: (In dj(k)) versus t for an initial state |(o;,5); 0) with 7/g = 5 and 
v = 5. The solid line corresponds to d em h = 6. The dotted lines correspond to 
values of d em i, from 7 to 10. 
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Figure 6.6: (In dj(k)) versus t for an initial state \a ; 0) with j/g = 5 and v = 10. 
The solid line corresponds to embedding dimension equal to 5. The dotted lines 
correspond to values of d em b from 6 to 10. 



cases. In contrast, there is no chaotic behavior for these values of the parameters 
if the initial state of the field is coherent. However, with a further increase in u, 
this situation changes, and (N) varies chaotically even for an initial coherent field 
mode. This is demonstrated in Fig. |T6] showing (lnd,(A;)) versus t for an initial 
state \a; 0) with v = 10. The corresponding value of A max is also quite large, 
« 0.60. 

To summarize: a time series and power spectrum analysis shows that, for 
small values of the ratio of the anharmonicity parameter to the strength of the 
inter-mode coupling, the dynamics of the expectation value of the mean photon 
number (or mean energy) of the field mode ranges from periodic to ergodic, but 
not chaotic, essentially independent of the nature of the initial state. On the 
other hand, for sufficiently large values of this ratio, the ergodicity properties of 
the sub-system representing the field mode depend significantly on the extent of 
departure from coherence of the initial field mode. In particular, A max is larger in 
the case of an initial field state that is a photon-added coherent state, compared 
to its value for initial coherent state. Table 6.1 provides a bird's eye view of these 
conclusions. 
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Increasing departure from coherence - 

Initial state 


a ; 0> 


(a, 1) ; 0> 


(a, 5) ; 0> 


Increasing nonlinearity 


-2 

Y/g=10 
|oc| 2 =v=l 


regular 


regular 


regular 


Y/g=10 
v=5 


regular 


regular 


regular 


y/g=5 

v=l 


regular 


regular 


chaotic 
*w -0.50 

max 


T/g=5 
v=5 


regular 


chaotic 
^max= " 58 


chaotic 


T/g=5 
v=10 


chaotic 


chaotic 


chaotic 

X ~ 1.00 

max 



Table 6.1: Qualitative dynamical behavior of the mean photon number of a single- 
mode electromagnetic field interacting with a nonlinear medium. "Regular" =>- 



A, 



0. 



We have already shown that the classical counterpart of the Hamiltonian under 
study (Eq. Q6.3Q ) corresponds to bounded motion in phase space, and that the 
system is integrable, with all its Lyapunov exponents equal to zero. We have also 
seen that the quantum Hamiltonian (Eq. (|6.1| )) commutes with the total number 
operator N to t , and has a discrete energy spectrum for all values of the parameters 
in H. What is the interpretation, then, of a positive value for A max as deduced 
from the time series for the sub-system variable (a^a) ( equivalent ly, (tfb))! 

This is best understood by appreciating the role of recurrence theorems and 
ergodicity in both classical and quantum dynamics. In classical physics, ergodicity 
itself is intimately linked with the Poincare recurrence theorem. This theorem es- 
sentially states that any phase space configuration of a system enclosed in a finite 
volume will be repeated to an arbitrary degree of accuracy after a finite, though 
perhaps very long, time. Quantum analogs of this classical recurrence theorem 
have been established and applied to diverse situations. Exact revivals or near- 
revivals of a quantum state occur [6|]] in systems with discrete energy eigenvalues. 
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Further, it has been argued |)5| that under any time-periodic Hamiltonian, a non- 
resonant, bounded quantum system will return to its original form several times 
during its dynamics, and this has been illustrated using computer experiments. 
Further, the link between wave packet revivals, and geometric phases in quantum 
physics and classical recurrences on the circle map has been examined p6| . In- 



spired by the seminal work of Fermi, Pasta and Ulam [67| on the recurrences of 



the initial state of a chain of nonlinearly coupled oscillators, several investigations 
have been carried out over the years on a wide variety of physical systems, in order 
to understand the conditions under which recurrences occur. In particular, it has 



been argued |68| that wave packet recurrences (exact revivals and near-revivals) 
and the Fermi- Past a- Ulam recurrences have the same origin: if the initial quan- 
tum state can be represented as a superposition of eigenstates, and the spectrum 
is non-equispaced, recurrences are bound to occur. 

On the other hand, high-resolution tracking of quantum and classical evolu- 
tions of a system will reveal differences in the dynamics of the expectation values 
of corresponding observables in these two cases. As indicated in Chapter 1, the 
origin of this dichotomy can be traced back to the inadequacy of the naive Ehren- 
fest theorem (which does not generally take into account the non-commutativity 



between x and p, as seen from Eq.(L6)) in retrieving the classical regime of the 
quantum system. An outcome of this feature is that the Lyapunov exponents that 
characterize the dynamical behavior of classical and quantum expectation values 
of the same observable can indeed be very different from each other |69 |. Two 



interesting observations [70] have been made in this context: (a) In isolated quan- 
tum systems with a discrete energy spectrum, using unitarity and the Schwarz 
inequality, it can be established that the Lyapunov exponents would vanish, if 
computed from time-series data collected over a sufficiently long time (which, in 
practice, could even turn out to be considerably longer than the characteristic 
time scales in the problem), indicative of non-chaotic behavior, (b) However, once 
measurement upon the system is included through appropriate interaction with 
an external system, the corresponding Lyapunov exponents need not vanish. This 
would suggest that, in principle, the dynamics of a subsystem, as deduced from 
time-series data of certain variables, may show exponential instability, even if the 
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system as a whole does not. The presence of hyperbolic (sets of) critical points 
in the corresponding classical version of the dynamical system plays a significant 
role in this regard. The "chaotic" behavior displayed by the mean photon number 
(a^a) in the system we have considered may be understood in the light of these 
observations. Ultimately the exponential instability associated with a positive 
Lyapunov exponent is an indication of the manner in which an initial wave packet 
spreads as it evolves in time. 
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CHAPTER 7 



Conclusion 

In the preceding chapters, we have described our results in detail, and also sum- 
marized them at appropriate places throughout the text. It remains to place the 
work in a broader perspective, and to list interesting open problems and avenues 
for future work. 

In general terms, the present work may be viewed as addressing an aspect of 
the broad theme of understanding the manner in which the notions of ergodicity 
and varying degrees of randomness, that are well-defined and notable features 
of classical dynamics, get modified in quantum dynamics. In specific terms, we 
have studied several aspects of wave packet dynamics that shed light on non- 
classical effects such as revivals and fractional revivals. As emphasized in the 
preceding chapters, such effects are striking manifestations of quantum mechanical 
interference phenomena. They show up in distinctive fashion in the expectation 
values of observables and their higher moments, in different kinds of squeezing, 
and in other indicators such as entropies and quasi-probability distributions, as 
elaborated upon in the foregoing. We have worked in the framework of two specific, 
rather simple, model Hamiltonians pertinent to a nonlinear medium. However, 
these models capture many of the essential features of nonlinearity, anharmonicity 
and inter-mode coupling. In this sense, the results arrived at enable us to draw 
conclusions that may be expected to remain valid in their broad features, if not 
in detail, in other modified situations as well. 

Our results also suggest a number of avenues for further exploration. We list 
some of these here, roughly in the order of progression of the text of this thesis. 

We have examined the revival phenomena exhibited by a single-mode coher- 
ent state or a photon-added coherent state as it propagates through a Kerr-like 
medium. It would be of interest to investigate their counterparts for other initial 



states, such as photon- subtracted coherent states and various squeezed states 
of light. 

The Wigner function is only one particular way of arriving at a non-classicality 
indicator in the situations considered in this thesis. Others, such as the Sudarshan- 



Glauber P-representation ]32] and the Husimi-Kano function [71], are also worth 
exploring, as is the extension of these considerations to the interacting two-mode 
case. 

The interacting two-mode case we have considered involves a bipartite Hamilto- 
nian. What would happen if this were not the case — for instance, if we considered 
three sub-systems coupled to each other? Further, if such a system were classically 
non-integrable, what sort of dynamical behavior would the quantum mechanical 
counterpart exhibit? 

The sub-system entropies we have considered (the von Neumann entropy and 
the linear entropy) are just two of several possible measures of entanglement, and 
they seem to exhibit essentially similar behavior, at least in the cases we have 
examined. What about other measures of entanglement, and what is the most 
appropriate or optimal measure for the purposes at hand? 

Finally, as perhaps the most interesting general question that arises in the 
present context: what are the possibilities and implications for quantum compu- 
tation using continuous variables (e. g., Gaussian wave packets), exploiting revival 
phenomena? It is evident that there exists a rich vein of interesting problems for 
further research in this area. 
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APPENDIX A 



Calculation of (a j(r a r+s } for an initial PACS |a,m) 
propagating in a Kerr-like medium 



The Hamiltonian governing the dynamics of the wave packet is given by H = 
%N(N — 1) where N = a^a. The initial state 1^(0)) is the PACS \a,m). The state 
at any t > is given by 



\i/>(t)) = e- l * N{N ~ 1)t \a,m). 



(A.l) 



Introducing complete sets of Fock states YlHo |/) (l\ (= 1) at appropriate places, 
we have 



oo oo oo oo 



1=0 n=0 u=0 v=0 

x (n| a) r a r+s \u) (u\ e -*xN(N-i)t |^ ^ ^ _ ( A2 ) 
Since our initial state is the PACS \a,m), it is convenient to use the notation 



{a^a^) m =m)\^a r+s m))- 



(A.3) 



When m = (corresponding to an initial CS), we will write the expectation value 
as simply (a) r a r+s ^j. On using 



a p \u) 



ill 



(u — p)\ 



1/2 



| It — p) 



and the orthogonality property of the basis states, Eq. (|A.2|) simplifies to 



(a^ r a r+s ) 
\ it 



1=0 



n 



V2 /! 



(l-r)\ 



-i X [(l+s)(l+s-l)-l(l-l)]t 



x (a, m\l) (I + s\a, m) . 



(A.4) 



Here, and in the rest of this Appendix, it is understood that contributions from 
terms of the form l/(— n)\, where n is a positive integer, vanish. (Note that the 
gamma function T(l — n) diverges for n — 1, 2, • • • .) Expanding \a, m) in the Fock 
basis, we find 

{a,m\l) = ^ =^7T^-y, , (A.5) 
y/m\L m (-u) {l-m)\' 

where v = \a\ 2 and L m is the Laguerre polynomial of order m. We use the 



definitions 72 



LM = f:( m ) tl f = ^/Vo (A.6) 



n=0 



and 



L M = V C + S ) = tEL*!^-*) (A.7) 



for the Laguerre and associated Laguerre polynomials. Equation (|A.4j ) reduces to 



CO 



,,s , ,l—m 



(a tr a r+s ) = - : : e ~<"e 

x /m m\L m {—v) ^ (l — r)\(l + s — m)\ (l—m)\ 



/! (l + s)\ a s v l „-, x , t .,- U i r -2; x ist 
But 



(1 + S)\ ,,l-m_,,-s d fl+. 



{i + s — my. dv m 

so that 

/vSp-fp-'Xsfs-l)' ™ „-2ixlst Jm . rl r ij l \ 

(a ]r a r+s ) = - — [ v r+s — — ) (A 10) 

\ lm u s m\L m (-u) (I - m)\ dv m \ dv r )' y ' ' 

This can be re-expressed as 

fr a B e -u e -i X s{s-i)t ™ / m \ ,d n u r+s \ d m ~ n+r ™ e - 2i x lst v l 

^ a ' m ~ v s m\L m (-is) ^\nj\ dv n J du m - n+r 2-* (l—m)\ 

(All) 

Simplifying, we find (after some straightforward algebra) 

n s p -v -ix{s-l+2m)st Jm . Jr . 

U r a r+S ) = \ — , , - (v r+s -?-[v m exvive- 2 ™*))). (A.12) 
x /m v s L m {-v\m\ dv m \ dv r L v 7 V v 7 
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This can be re-cast in the form 



(a^a r+s ) = a e e y (>\ 



v s L m (— v) ^— ' \n J (m — r + n)\ 

n=n min 

x — [u m+s+n exp (ise- 2isx % (A.13) 



where n min = max (0, r — m). Using the definition in Eq. ( A.7 ) for the associated 
Laguerre polynomial L s m we get, finally, 

<a tr a r+s ) m = a s exp [ - v - t X {s -l + 2m)st + u e~ 2isxt ] 



r 



( r ) m: ^ c zi ^™ 11 . (A.i4) 



n) (m — r + n)\ L m (—u) 



This form of the expression for (a^ r a r+s ) m at an arbitrary time t, corresponding 
to the initial state \a,m), is used for numerical calculations. 



Setting m = in Eq. ( A. 14 ), we obtain 



( a t r a r+s ) = a s u r e- v ^- coa 2sxt) exp [-%x(s{s - 1) + 2rs) t - iv sin 2sx*] (A. 15) 
as the corresponding expression for an initial coherent state |a). 
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APPENDIX B 



Analytic expressions for (x(t)) and (p(t)) for an 

initial PACS 



Setting r = and s = 1 in Eq. (|A.14|) , we get 



(a) m = — e~ 2m ^ exp (ue- 2 ^ 1 ) L x m {-v e~ 2 ^). (B.l) 

The expectation values (x(t)) and (p(t)) m (i.e., the mean values of x and p at 
time t in the case of the initial PACS \a,m)) can now be calculated using Eq. 
( |B.1|) and its complex conjugate, which yields (at) . 

Recall that 



(x(t)) m = V2Re (a) m and {p{t)) m = V2 Im (a) m , (B.2) 



where we emphasize that (a) m is computed at time t, and the suffix m is meant 
to denote the fact that the initial state considered is \a,m). We find 



/2 e -i/(l-cos2x<) r . 

(z(t)) m = T , s Re aexpH(2m X t + ^sin2xt)]^(-z/e- 2 ^) 

(B.3) 

and 

/2 e -v(l-cos2xt) f . 

(P(t)) m = n S Im aex P H(2mxt + z/sin2 X t)]^(-^e- 2 ^) . 

(B.4) 



Setting m = in the above equations, we obtain the corresponding expressions 
for the case of an initial CS, as written down in Eqs. ( |2.19| ) and ( |2.20| ). 



Returning to the case of an initial PACS \a,m), we see from Eqs. flB.3|) and 



4]) that the initial values of the expectation values concerned are given by 

<*(0)> m = xo , <p(0)) m = po , (B.5) 

where x$ + zpo — Oi\/2. The difference 

(^)> m ^W)) m = ^« m (B.6) 
deserves special mention. We find 

«t)) m - * (P(t)) m = V2-^—-e 2 ^ exp {ue 2 ^) L^-ve**). (B.7) 
This can be rewritten in the compact form 

<s(*)L - < <p(*)L = (*o - ipo) ~ j { f U , (B.8) 

where ((t) = exp (2ixt). 

It is convenient to define the quantities 

X m (t) = (x(t)) rn exp [v (1 - cos2xt)] and P m (t) = (p(t)) m exp [v (1 - cos2x*)]. 

(B.9) 

Using the expressions for (x(t)) m and (p(t)) m in Eqs. ( |B.3| ) and ( |B.4j ), we get 

X m (t) = — ^ — -Re\( y x + ip )exp[-i( y 2mxt + usm2xt)}L 1 m (-ue~ 2ixt )}, 
P m (t) = 1 Im {(x + i Po ) exp [-i{2m X t + z/sin2 X t)] ^(-z/e" 2 ^)). (B.10) 



These can be written after some straightforward algebra in the concise form 

X m (t) = x Re z m (t)+ po Im z m (t) , 
P m (t) = Po Re z m (t) - x Im z m (£) , 
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where 



-v e 



2i x t\ 



Lmi—lA 



exp [i (2mxt + v sin 2%t)] 



(B.12) 



For ready reference, we have written out Eqs. (|B.11|) and (|B.12|) once again in 
Chapter 3 (Eqs. ( p.8|) and (|3.9|), respectively). In deriving the expressions above, 
we have made frequent use of the fact that e z and L s m (z) are real analytic functions 
of z: that is, they satisfy f(z*) = f*(z). 
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APPENDIX C 



Calculation of D q m (t) for an initial PACS \a,m) 



Setting r = and s = q (a non-negative integer) in Eq. (|A.14j) of Appendix A 
and using the series expansion for the associated Laguerre polynomial, we get 



(A 



cos 2q\t) 

L m {-v) 



n=0 



[m + q)\ v n 



m — n) \ (n + q) \ n\ 



x 



exp { — i [(q — 1 + 2m + 2n) xqt + vsm2xqt — qO]}, (C.V 



where 9 is the argument of a, that is, a = v x l 2 e %e . Similarly, when r = q and 
s = 0, Eq. (|A.14|) gives 



U q a q ) = — 

\ I m T 

Lit 



1 " 

f-iA 



ml 



n I [m — q + n)\ 



v n Ll(-v 



m\ 11 



(C.2) 



where n m [ n is now given by max(0,g — m). The quantity Dq(t) for an initial 
state | a, m) can now be calculated using the definition (see Chapter 4, Eq. ( |4.12|) ) 



Re{a 2q ) m -2(Re(a q ) m ) 2 + (a^a q ) j 



(C.3) 



Here -F 9 (N) = [a q , a^]. After some algebra, we find that the operator F q (U) is a 
polynomial of order (q — 1) in the number operator N = a^a, given by 



F q (H) = q\ 



n=l ^ ' 



(C.4) 



Since N commutes with the Hamiltonian H = \ N(N — 1), it follows that (F q {N)) 
is also a constant of the motion, and its value remains the same as its initial value 
computed in the state \a,m). 



Substituting the values for Re {a q ) m , Re (a 2q ) m and (a^a 9 )^ in Eq. ( |C.3| ), we 



get 



\ Lm (-u) (F,(N)>„ ,D<™ >(t) = f](™+£)^ 

n=0 



X COS 



2(2m + 2n + 2g - l)xgt + 1/ sin Axqt - 2q9 



2( ■-2"fl-n.s2\<//) ^ /// • 7 \ 



+ E f g V ^ r^ L^-v), (C.5) 

71 — 71 m j n 

where n min = Max (0 , q — m). When m = (corresponding to an initial CS), 
simplifies to 

D « ){t) = ^^{ 1 + e ~ 2 ^ 229Xtcos ( 2 XQ(^-^t + ^smA X qt-2 q e) 
- 2e- 4iysin2,?xt cos 2 (xq{q-l)t + usn2xqt-q6)y 

(C.6) 



As mentioned in the text, Dq°\t) has been written as simply D q (t) in Sec. 4.2, in 
the case of an initial CS. 



100 



APPENDIX D 



Density matrix for a field mode interacting with 

a nonlinear medium 



The Hamiltonian which describes the interaction of a single-mode field interacting 
with the nonlinear medium is given by 



H = Lua^a + Lu b j b + -ftf 2 b 2 +g(a f b + tfa). (D.l) 

The total number operator N tot = (a^a + tfb) commutes with H. We use N to 
denote its eigenvalues. The eigenstates \ipNs) of H are then labeled by integers iV 
and s, where N = 0, 1, ... ad inf. and s — 0, 1, . . . , N for a given value of N, 
as in Chapter 5, Sec. 5.2. These can be expanded in terms of the basis states 
\N — n) a ® \n) b = \N — n; n) as 

N 

\^Ns)=J2 d n S \N-n;n), (D.2) 

n=0 

with expansion coefficients 

d^ = (JV-n;n|^.). (D.3) 



The Hamiltonian in Eq. ( p.l|) can be represented in the basis \N — n\ n) as a 



real, symmetric, tridiagonal matrix. The eigenvalues Aat s and eigenstates I^tvs) of 



H can be found numerically using appropriate matrix algebra routines Jf3 
The unitary time evolution operator is given by 



oo N 



U(t) = Yl ^ (~ iX Nst) \*pNs) {lpNs\ ■ (D.4) 



N=0 s=0 



Hence, an initial state \ip(0)) evolves in time to the state 

oo N 

m)) = u(t) wo)) = e E ex p h-w) (i>m\m) i« . (d.5) 

7V=0 s=0 

The corresponding time-dependent density matrix is 



oo N oo AT' 



p(t) = EE E E exp [-z(Ajv s -AArv)^] (^1^(0)) (^(0)|^'s') |VW (V'JV's' I • 
N=0 s=0 N'=0s'=0 

(D.6) 

If the atomic oscillator is initially in the ground state |0) 6 and the field is in 
the Fock state \N) a , we have \if>(0)) = \N; 0), and 

N N 

p(t) = E E exp H( A ^ - *Ns>)t] (^ns\N; 0) (N ; 0|V>av) \^ Ns ) (V>av| ■ ( D - 7 ) 

s=0 s '=0 

(This expression is explicitly iV-dependent, as expected.) Substituting for d^ s 
from Eq. ( p.3|) , p(t) can be re-expressed in the form 

N N 

p(t) = E E ex p - A *«')*] d ° s C i^-> ■ ( D - 8 ) 



s=0 s'=0 



Hence 

(ip M i\p(t) \ipM'i>) = exp [-i(A M ; - A M 'Z') * ] <C rf o f ' r <W<W- (D.9) 

It is evident that is an (AT + l)-dimensional diagonal matrix in this case. 
For an initial state \(a,m)] 0) we find, using Eq. (p.6|) , 



°° N °° N ' VWW^. (a) N ~ m {a*) N '- m 



*m = mn e ( -,a E E E E 



m\L m (-u) ^ ^ ^ ^ (N-m)\ (N' - m)\ 

mv N=m s=0 N'=m s'=0 V 1 V ' 



x exp [-i(X Ns - \ N r s/ )t)d$ s d£' s ' \ip Ns ) (V>jvv| , (D.10) 
where we have used Eq. Q3.1]) of Chapter 3, and expanded |a, m) in the Fock basis. 
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The corresponding matrix elements of the density matrix are given by 

e~ v VM\ M'\ (a) M - m (a*) M '- m 



(4>Ml\p(t) \lpM<l<) 



m\L m {-v) (M — m)\ (M' — m)\ 
x exp [-i(X Ml - \ M , v )t)dtf l df v . (D.li; 



Here, and in the rest of this Appendix, it is understood that contributions from 
terms of the form l/(— n)\, where n is a positive integer, vanish. 

The expectation values and higher moments of the quadrature variables £(£) 
and r)(t), defined in Eq. (|5.9|) of Chapter 5, can now be obtained numerically, 
using the above expressions for the density matrix and the matrix elements of the 
operators a and b in the basis \ipNs)- The latter are given by 

N> 

ftM a |Vivv> = ^2(N' - nf' 2 < s d%' s ' 5 N , N ,-i (D.12) 

n=0 

and 

N' 

(lpNs\ b |^W> = Yl 721/2 d n-l d n S> 5 N,N'-1 , (D.13) 
n=X 

respectively. As these are purely off-diagonal, and p(t) is diagonal for an initial 
state that is a direct product of Fock states, it follows that all the odd moments 
of £ and 77 vanish identically for all t, as asserted in the text. This is no longer 
true for the other classes of initial states considered. 

The time-dependent reduced density matrices pkif) (k = a,b) are given by 
p a (t) = Tr b \p(t)] = ^2 b (n\p(t)\n) b , 

n=0 

00 

p b (t) = Tr a [p(t)] = J2a(n\p(t)\n) a . (D.14) 

n=0 

Corresponding to an initial state \N ; 0), these reduced density matrices pk(t) take 
the form 

N N N 

Pa(t) = exp[-i(\ Ns - \ Ns >)t] 

n=0 s=0 s'=0 

x d» s C d N J \{N-n)) a a {(N-n)\ (D.15) 
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and 



N N N 



n=0 s=0 s'=0 

x d^d^' d^_ n d^: n \(N-n)) bb ((N-n)\. (D.16) 



Hence we have, in the Fock basis, 

N N 

a(n\p a {t)\n') a = ^^eX P H(A^-A^0i]CC'<in4-n^nn' 
s=0 s'=0 

(D.17) 



and 



N N 

M Pb(t) \n) b = £ £ exp H( A "» " A ^') * ] d o S <C C d% s ' 5 nn , . 

s=0 s'=0 

(D.18) 



As before, these are explicitly iV-dependent finite-dimensional matrices. 
For an initial state \(a,m) ; 0) the expressions for puif) are given by 

e~ v ^ A A ^ v/aTTM (a) N ~ m (a*) N ' 



m\L m (-v)^ A^ A^ A^ A^ (N-m)\(N'-m)\ 

mv I n =0 N=N min s=0 N'=N min s'=0 V 1 V 7 

x exp [-i(Ajr. - \ N , S ,) t ] C < v C df' K^V - n)) a a ((JV' - n)| 

(D.19) 



and 



oo oo N oo AT' 



^ = rn ir...r-,^ 2^ 2. Z, Z, 



m!L ro (-z/)^ A^ A^ A^ A^ (N-m)\(N'-m)\ 
x exp [-i(Ajv s - \N>s')t]d% s dp' d^ n d% s L n \(N - n)) b b ((N' - n)\ , 

(D.20) 

where iV m i n = max (n,m). These are infinite-dimensional matrices. The corre- 
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sponding matrix elements of pk{t) in the Fock basis are given by 

e -» y ^ ^ V(n + 0!(n + Z')! («) n+ '" m («T +r 



m\L m {-v) ^j^f^ (n + I - m)\ (n + I' - m)\ 

x exp [-z(A (n+0s - A ( „ w ) * ] ^ 4" +0s 4 n+0s ' 

(D.21) 



and 

e- yy y >/(n + 0!(n + IO! (a)^ (q*)"+''-" 



(n+l)s ,(n+l')s' ,(n+l)s A (n+l')s' 



x exp [— i(\n s — Ajv's') t] d d d\ d\, 



(D.22) 



where n min = max (m — l,m — I'). 



All the reduced density matrices (with elements given by Eqs. ( 1LH7I) , ( |DjJ) 



( P-21|) and ( p.22|) ) are hermitian. They are diagonalized numerically, and their 



eigenvalues are used to compute the entropies Skit) and Skit) given by Eqs. Q5.8|). 
In the case of infinite-dimensional matrices, convergence in numerical computation 
is provided by the factorials in the denominators of the summands in the expres- 
sions derived above for the matrix elements. We use double precision arithmetic 
with an accuracy of 1 part in 10 6 . As mentioned in the text, we use the equality of 
p a {t) and pb{t) as one of the checks on the numerical computations. Other checks 
include the condition p 2 (t) = p(t) for the density matrix of the total system, as 
we are only dealing with pure states. 
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